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ABSTRACT 
In  the  modern,  high  performance,  supersonic  aircraft,  varying 
conditions  of  flight  produce  airframe  transfer  functions  whose  para- 
meters' also  vary.  Fixed  parameter,  servo  control  loops  in  general  can- 
not provide  satisfactory  control  over  the  entire  region  of  flight.  Self- 
adapting  control  loops  which  automatically  alter  their  parameters  in 
order  to  provide  the  most  satisfactory  control  response  are  becoming        y 
increasingly  important.  The  analysis  of  several  lateral  airframe 
responses  in  order  to  determine  which  is  most  favorable  for  adaptive  con- 
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trol  is  discussed.  The  design,  analysis  and  description  of  a  method  for 

II 
self-adaptive  control  of  the  yaw  augmentation  control  loop  of  a  high  per- 
formance aircraft  are  then  discussed  with  appropriate  conclusions  antf  . 
recommendations  as  to  further  research  or  development. 
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Chapter  I 
Introduction 

In  the  development  of  modern  weapons  systems ,  during  the  years 
following  World  War  II,  certain  rapid  strides  in  engineering  techniques 
as  well  as  technologies  have  resulted  in  weapons  concepts  which  were  here- 
tofore unimaginable.  The  development  of  the  jet  engine  for  aircraft  and 
the  rocket  engine  for  missiles  or  experimental  aircraft  have  created 
demands  for  system  control  in  the  regions  and  environments  of  flight 
which  occur  at  velocities  several  times  the  speed  of  sound. 

With  the  advent  of  the  supersonic  aircraft  and  its  increased  opera- 
tional capabilities,  numerous  problems  never  before  encountered  in  such 
magnitude  and  importance  are  foremost  in  the  designer's  «ind.  He  must 
reckon  with  aircraft  velocities  which  are  subsonic  as  well  as  supersonic, 
and  dynamic  pressures ,  q ,  whose  magnitudes  range  from  those  normally 
encountered  with  conventional  aircraft  to  those  associated  with  smaller 
types  of  missiles  and  rockets .  The   designer  must  take  into  account  the 
high  maneuvering  rates  possible  under  such  flight  conditions  and  the  in- 
herent sensitivity  of  such  an  aircraft  to  control  surface  forces  in  the 
extreme  conditions  of  flight.  Control  measures  must  be  supplied  which 
will  operate  to  minimize  the  large  forces  which  could  occur  on  the  air- 
crafts'  structure  and  cause  its  destruction.  In  addition,  these  control 
measures  must  account  for  the  human  factors  which  unfortunately  enter  the 
picture  in  high  performance,  high  velocity  situations.  Human  response  time 
to  instantaneous  accelerations  and  disturbance  forces  is  prohibitively 
slow.  The  human  is  also  relatively  insensitive  as  to  the  proper  magnitude 
of  an  input  command  to  the  control  surfaces.  In  a  high  performance  super- 
sonic aircraft,  these  improper  inputs  could  result  in  large  forces  which 
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might  exceed  the  structural  load  limits,,  In  addition,  the  lag  in  human 
response  time  could  result  in  corrective  action  reinforcing  a  rapid  oscil- 
latory disturbance  so  as  to  create  less  stable,  or  perhaps  unstable,  oper- 
ation. 

The  inherent  increase  in  the  performance  capabilities  briefly  mentioned 
has  created  additional  problems  in  terms  of  aircraft  stability  and  response . 
No  longer  does  the  aircraft's  dynamic  response  vary  slightly  under  all 
possible  conditions  of  flight.  The  variance  in  dynamic  response  at  ex- 
treme conditions  of  flight,  i.e.;  landing  characteristics  as  compared  to 
high  altitude  high  velocity  flight,  is  large.  The  aircraft  may  conceivably 
exhibit  good  stability  undergone  flight  environment  and  at  another  altitude 
and  velocity  the  airframe  may  become  unstable.  Clearly,  such  an  aircraft 
dominates  the  pilot's  capabilities  to  control  and  command.  There  is 
established  a  need  for  some  type  of  automatic  control  to  assist  the  pilot 
in  his  function. 

Automatic  control  systems  were  developed  to  provide  stable  aircraft 
response  characteristics  in  the  pitch,  yaw  and  roll  modes  of  motion.  The 
basic  requirements  being  that  the  autopilot  provide  the  ability  for  the 
aircraft  to  remain  steadily  on  the  commanded  course  or  maneuver,  while 
providing  the  necessary  stability  in  the  response  modes  mentioned  above  0 
These  requirements  pertain  to  any  disturbances  which  may  occur  as  inputs 
to  the  airframe .  Such  inputs  may  include  undesirable  oscillations  or 
rates  due  to  commanded  inputs ,  gust  disturbances  and  random  force  inputs , 
such  as  buffeting,  in  the  transonic  region. 

Generally,  autopilots  augment  the  commands  generated  by  the  pilot. 
The  pilot  possesses  complete  control  of  the  aircraft's  command  inputs. 
The  autopilot  in  this i instance -serves  to  provide  stable  airframe  maneuver 
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to  the  desired  attitude  and  heading.  Recently,  however,  in  the  develop- 
ment of  automatic  control  systems  for  supersonic,  high  performance  air- 
craft, certain  command  functions  normally  handled  by  the  pilot  have 
necessarily  been  programmed  into  the  autopilot  which  automatically 
applies  the  aircraft  command  to  achieve  the  desired  response.  This 
paper  concerns  itself  with  such  an  instance  in  the  yaw  control  channel 
of  the  autopilot  for  supersonic  aircraft.  ^^ 

As  previously  indicated,  the  aerodynamic  characteristics  of  this 
particular  type  of  aircraft  vary  over  the  entire  region  of  flight  capa- 
bilities. There  must  be  some  compromise  in  the  designer's  planning  to 
account  for  the  necessary  stability  required  at  landing  as  well  as  at 
high  velocity,  high  altitude  flight.  Obviously  the  physical  dimensions 
of  the  aircraft  cannot  be  modified  so  as  to  always  present  an  ideal  design 
with  respect  to  air  flow  at  any  and  all  flight  conditions.  It  follows, 
then,  that  the  automatic  control  system  must  be  capable  of  providing  a 
control  which  is  satisfactory  for  all  of  the  flight  conditions ,  and 
therefore  any  changing  aircraft  characteristics  which  may  occur 0  In  many 
instances  this  control  is  satisfactorily  achieved  by  servo  loops,  con- 
taining constant  parameters ,  which  provide  acceptable  response  regardless 
of  the  flight  condition.  However,  in  the  case  of  the  aircraft  considered 
by  this  report,  a  constant  parameter  control  loop  was  not  sufficient  to 
insure  acceptable  response  throughout  the  regions  of  flight.  A.  variation 
of  one  or  more  parameters  of  the  control  loops  were  considered  necessary 
to  provide  the  desired  stability  requirements.  Such  a  control  loop  would 
then  adopt  itself  to  each  flight  condition  by  sensing  its  inadequacy  and 
providing  the  necessary  changes  in  its  parameters  to  improve  the  response. 

£  J      Shall  refer  to  references  listed  at  the  end  of  this  report. 


The  research  on  the  adaptive  controller,  which  refers  to  this 
entire  unit  with  the  capability  of  varying  control  system  characteristics, 
has  been  carried  out  principally  in  recent  years 0  The  United  States  Air 
Force  became  interested  in  the  development  of  advanced  aircraft  design 
which  indicated  a  need  for  relatively  more  sophisticated  control  methods. 
Soon  leading  laboratories  in  the  control  systems  development  field  were 
contributing  to  the  problems'  solution.  Many  approaches  were  made  to  the 
problem.  Some  strived  to  obtain  satisfactory  response  by  single  parameter 
variation  as  a  function  of  perhaps  mach  number,  M,  and  altitude,  h.  In 
such  a  method  some  switching  was  incorporated  to  set  a  predetermined  loop 
gain  which  w  ould  provide  satisfactory  airframe  response  over  some  region 
defined  by  limits  expressed  as  some  function  of  M  and  h.  Other  systems, 
certainly  more  sophisticated,  used  the  model  reference  system  whereby  an 
analog  model  of  the  aircraft  received  the  same  inputs  as  the  airframe. 
By  measuring  the  airframe's  response  and  comparing  this  to  the  model's 
response,  which  was  considered  the  ideal  response,  an  error  signal  was 
generated  which  served  to  provide  some  measure  of  the  system  changes 
necessary  to  match  the  physical  to  the  ideal  model  response.  This  error 
signal  was  handled  in  numerous  ways  by  the  research  and  development  units . 
By  developing  and  applying  criteria  of  optimum  response  to  this  error, 
input  signals  could  be  processed  and  shaped  so  as  to  provide  forcing 
functions  whi  ch ,  when  applied  to  the  airframe ,  caused  it  to  reproduce  the 
model's  response;  or,  the  error  criterion  could  provide  the  impetus  for 
varying  loop  gains  and/or  time  constants .  Clearly  there  are  many  varia- 
tions and  applications  of  these  and  other  principles  not  discussed.  From 
the  summary  of  these  systems ,  it  can  be  readily  appreciated  that  there 
are  requirements  for  accurate  data  input  in  order  to  predict  switching 
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function  values,  usually  provided  by  air  data  computers,  or,  analog  com- 
puters are  necessary  to  provide  the  model  reference  and  integrate,  in 
cases,  functions  which  modify  the  input  signal1 s  shape.  Digital  computers, 
by  numerical  methods,  could  handle  some  of  these  analog  functions  but  the 
need  for  some  computer  remains. 

Many  of  these  systems  have  been  satisfactorily  flight  tested  or  under- 
gone extensive  tests  under  simulation  conditions 0  Applications  of  adaptive 
controls  in  processes  are  becoming  prevalent  in  industry  as  automation 

becomes  more  evident.  Their  future  application  in  fields  where  human 

^2,3,4,7,8,97 
inadequacies  are  detrimental  to  control  is  insured. 

The  material  presented  by  the  authors  in  this  thesis  is  the  analysis 
of  a  high  performance  and  velocity  aircraft  with  intentions  of  providing 
a  satisfactory  response  over  all  flight  conditions  by  simply  the  intro- 
duction of  a  gain  change  in  the  control  loop.  Furthermore,  the  impetus 
for  furnishing  the  gain  change  is  to  be  entirely  derived  from  the  airframe 
response  and  is  therefore  a  self  adapting  scheme.  No  exterior  inputs,  such 
as  M  and  h,  were  considered.  No  computers  were  utilized  in  the  scheme. 
The  aircraft  response  was  to  approach  the  "selected  ideal"  as  closely  as 
possible  in  order  to  define  the  optimum  response. 

The  study  was  carried  out  at  the  Autonetics  Division  of  North  American 
Aviation,  Incorporated  under  the  direction  of  Mr.  Richard  K.  Smyth,  senior 
research  engineer  and  supervisor  of  the  Advanced  Analysis  Unit  of  the 
Flight  Control  Section.  The  study  was  conducted  in  June  and  July,  1959. 


Chapter  II 
General  Approaches  to  the  Problem 

The  development  of  the  adaptive  controller  as  an  integral  part  of  an 
automatic  flight  control  system  has  been  briefly  mentioned  in  Chapter  I. 
Generally  the  systems  evolved  from  basically  simple  ideas  and  techniques 
associated  with  obtaining  the  "optimum"  performance  from  a  servomechanism 
or  control  system.  Just  what  is  meant  by  the  term  " optimum"  depends 
primarily  on  the  user  of  the  word.  Reference  5  uses  the  word  optimum  to 
define  the  most  commonly  acceptable  second  order  system  damping  ratio, 
0.7.  Still  other  authors  define  the  system  as  optimum  when  certain  math- 
ematical relationships  are  introduced  which  alter  the  forcing  function, 
and/or  the  system  parameters,  giving  a  "desired"  response.  An  example 
of  this  would  be  the  least  squares  criterion  proposed  in  reference  6. 
It  seems  evident  then  that  the  field  of  adaptive  servome onanisms  and 
optimum  design  research  are  interwelated .  Furthermore,  the  basis  of 
adaptive  control  is  some  method  of  attaining  the  most  desirable  response 
under  variable  operating  conditions  or  inputs „  Hence ,  the  most  desirable 
response  is  often  referred  to  as  the  "optimum  response"  under  these 
imposed  conditions.  Therefore,  it  is  a  recognized  fact  that  the  w state 
of  the  art"  of  adaptive  mechanisms  to  date  have  relied  upon  these  techniques. 

The  first  portion  of  this  report  concerns  itself  with  the  analysis 
of  several  recommended  airframe  lateral  responses  with  the  object  in  mind 
that  these  chosen  transfer  functions  satisfy  the  following  requiremstita . 
First,  the  response  must  furnish  oscillatory  stability  over  the  entire 
region  of  flight.  Furthermore,  any  divergent  instability  from  a  real  root 
must  be  slight  enough  to  permit  the  pilot  or  autopilot  to  remain  on  course 
by  command  inputs.  Secondly,  the  response  characteristics  should  provide 
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some  similarity  to  the  response  of  a  second  order  system  in  order  that  the 
human  pilot  preference,  3  rad/sec  and  a  damping  ration  of  0.7,  as  well  as 
second  order  optimizing  principles  might  be  attempted.  ^   '  "* 

The  recommended  responses  to  be  investigated  as  indicative  of  lateral 

aircraft  motion  were  lateral  acceleration  of  the  aircraft  center  of  gravity, 

d*P  d*  0     • 

IL,  yaw  rate,  ,f      or  r,  and  yaw  acceleration,  dJ^  or  r,  in  combination 

with  r  or  N  .  These  functions  are  the  descriptive  characteristics  of  an 
airframe's  lateral  motion  and  are  capable  of  measurement  by  instrument- 
ation perfected  and  existing  as  useable  hardware.  Pig.  1  indicates  the 
basis  for  measurement  of  the  quantities  by  relation  to  the  basic  stability 
axes  contained  in  the  airframe .with  their  origin  at  the  e.g.  or  center  of 
gravity.  Positive  N  is  measured  along  the  positive  y  axis  and  is  the 
result  of  a  left  rudder  deflection 0  Positive  y    ,  r  and  r  are  measured 
in  a  clockwise  direction  about  the  vertical  axis . 

The  stability  analysis  was  carried  oat  with  the  following  basic 
assumptions .  The  analysis  was  to  be  linear  with  all  responses  as  a  function 
of  rudder  deflection  only.  All  other  control  surfaces  were  assumed  to  re- 
main at  zero  deflection.  The  aircraft  was  allowed  three  degree  freedom  as 
a  result  of  these  rudder  deflections  and  was  assumed  initially  in  steady 
level  flight  prior  to  application  of  a  forcing  function.  Six  indicative 
flight  conditions  representing  the  extremes  of  the  designed  region  of 
flight  were  chosen  as  the  basis  of  requirements  for  stability.  These  six 
flight  conditions  represent  individual  transfer  functions  and  therefore 
varying  responses. 

The  basic  control  loop  is  shown  in  Fig.  2.  The  loop  consists  of  the 

transfer  functions  for  the  airframe,  necessary  instrumentation,  compensation 

if  required,  and  provision  for  variable  gain  adjustment. 

To  summarize,  the  analysis  was  necessary  in  order  to  determine  which 
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of  the  lateral  aircraft  responses  at  extreme  conditions  of  flight  provided 
a  response  to  rudder  deflection,  only,  which  was  stable  in  oscillation  and 
approximated  a  second  order  system  at  all  six  flight  conditions.  This 
was  necessary  to  permit  application  of  basic  optimum  criterion  to  provide 
satisfactory  adaptive  control. 

The  design  of  the  adaptive  control  circuits  was  based  on  the  fol- 
lowing desired  characteristics.  The  designed  scheme  was  to  be  simple  and 
easily  mechanized.  The  scheme  was  required  to  base  its  own  adjustments 
on  an  error  criterion  which  exhibited  selectivity  of  an  acceptable  damping 
ratio  based  on  information  obtained  from  the  control  loop  and  not  from 
sources  external  to  the  loop.  The  adaptive  controller  was  to  achieve 
this  damping  ratio  by  a  gain  variation  only.  Therefore,  effectively  the 
adaptive  control  system  would  be  capable  of  sensing  the  requirements  for 
stable  response  by  looking  at  the  contemporary  airframe  dynamics. 

The  various  chapters  contained  within  the  Scope  of  This  Investigation 
are  logically  arranged  to  provide  introductions  and  explanations  of  the 
sequential  development  of  the  stability  analysis,  self  adaptive  design  and 
test  of  the  entire  system.  Specific  Appendices  are  referenced  as  their 
need  arises  in  the  report  text.  Conclusions  of  this  report  and  reconnend- 
ations  for  further  design,  test  and  analysis  are  presented  as  a  final  summary 
of  the  material. 
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Chapter  III 
Scope  of  This  Investigation!  Stability  Loop 

Referring  to  Pig.  2,  it  is  considered  appropriate  to  detail  the  make- 
up and  function  of  the  loop  parameters.  The  basic  control  loop,  as  briefly 
explained  in  Chapter  II,  serves  as  the  yaw  augmentation  loop  for  the  air- 
craft^ autopilot  or  automatic  flight  control  system.  There  are  no  external 
command  inputs  to  the  yaw  loop  for  this  purpose .  Therefore ,  the  loop 
becomes  one  in  which  the  airframe  is  disturbed  by  some  force  which  imparts 
a  measurable  lateral  response.  The  measured  response  is  fed  back,  with  the 
appropriate  sign  of  feedback,  to  actuate  the  rudder  in  order  to  reduce  the 
response  to  zero.  Since  the  control  loop  is  a  single  closed  loop,  and 
referring  to  Fig.  3,  we  may  combine  the  loop  transfer  functions  into  a 
single  function  with  a  variable  gain  factor  and  unity  feedback.  Pig.  U 
indicates,  however,  that  to  obtain  the  response  of  the  airframe,  as  an 
example,  the  problem  reverts  to  one  in  which  there  is  a  direct,  feed  for- 
ward function  and  a  feedback  function. 

The  instrumentation  contained  within  the  control  loop  was  as  follows; 

a.  autopilot  actuator 

b.  main  rudder  actuator 

o.  measuring  instrument  (rate  gyro,  accelerometer ,  or  angular 
acceleromster) 

The  autopilot  actuator  was  a  closed  loop  servo,  hydraulic  control  valve 

which  regulated  the  hydraulic  pressure  fed  to  the  main  rudder  actuator. 

The  main  rudder  actuator  controlled  the  displacement  of  this  vertical 

control  surface.  The  servo  valve  used  in  the  autopilot  actuator  exhibited 

a  closed  loop  second  order  response  function  of 

°c<8>        .  800 

OrCb)  AA         (s  ♦  20  t  J20 
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while  the  main  rudder  actuator  had  a  response  characterized  by  the  function 

g^TsT  m  "  (s  +10) 

Investigation  of  these  functions  by  an  application  of  the  final  value 

theorem,  for  any  particular  input,  shows  that  for  the  autopilot  actuator 

where  the  input  function  is  given  by  QR(s)  and  the  output  function  by 

©  (s),  we  find  that 
c 

0„  800VB>AA         _  8008r(8>AA 


c       AA       (s  ♦  20  ±  J20)  s*  ♦  4-Os  ♦  800 


then; 


800  0-(s) 
lim  9c(t)  =  lim  s.Oc(s)  =  lim  s  •  g>  ,  y*   ,  800 

t  -*oO       s  ->0      S  -+  O 

»  lim    s  •  OoCs) 
s-*0      R 

and  observe  that  the  steady  state  value  of  the  instrument  is  unity.  This 

also  holds  true  for  the  main  rudder  actuator  and  the  measuring  instruments . 

The   measuring  instruments  used  in  the  analysis  were  a  rate  gyro  in  order  to 

measure  yaw  ratef  an  accelerometer  to  measure  lateral  acceleration  at  the 

c0g.  of  the  airframej  and  an  eccentric  angular  accelerometer  developed  by 

North  American  Aviation  Corporation  which  was  capable  of  measuring  the 

lateral  acceleration  at  the  e.g.  and  some  desired  proportion  of  r.  The 

signals  were  essentially  combined  to  give  a  single  output  value  from  the 

device .  In  each  instance  the  identical  transfer  function  was  used  since 

it  provided  a  measure  of  simplicity  and  was  a  fair  approximation  to  all 

three  instances.  The  transfer  function  used  was 

Ms) 


-        15795 
"o^TsT     (s  +87.9  i  J89.7) 
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and  is  second  order  with  a  steady  state  value  of  unity.  Its  dynamic 
characteristics  were  a  natural  frequency  of  20  cps  and  a  damping  ratio 
of  0.7. 

In  addition  to  the  instrumentation  contained  in  the  loop  was  the 
appropriate  transfer  function  of  the  airframe  as  a  function  of  rudder 
deflection,  b_.  A  particular  aircraft  function  was  inserted  into  the  loop 
equation  for  each  of  the  desired  lateral  responses  at  every  condition  of 
flight.  This  required  analysis  then  of  a  minimum  of  eighteen  airframe 
functions.  These  functions  will  be  given  and  explained  in  detail  in  the 
following  Chapters. 

Finally,  the  loop  provided  for  the  insertion  of  recommended  compen- 
sation, a  variable  gain  adjustment  and  some  means  for  reversing  the  phase 
of  the  loop  signal,  if  required 0 


15 


Chapter  IV 
Scope  of  this  Investigation:  Lateral  Acceleration 
as  a  Dynamic  Response 

By  referring  to  Fig.  2  which  indicates  the  control  loop  for  yaw 
augmentation  the  transfer  functions  for  lateral  acceleration  as  a 
function  of  rudder  deflection,  Hj/&d  »  ma7   be  inserted  as  the  particular 
aircraft  function.  Since  it  is  customary  to  express  the  function  £u  in 
degrees,  the  gain  variable  in  the  loop  will  be  expressed  as  degrees/g 
and  designated  K-  . 

The  transfer  functions  of  N/6u  a*  the  specific  conditions  of 
flight  are  listed  in  Table  T,     The  extreme  flight  conditions  chosen  to 
represent  the  limits  of  aircraft  performance  are: 

Condition  1:  low  altitude,  low  dynamic  pressure,  landing  condition 

Condition  2%     low  altitude,  low  dynamic  pressure 

Condition  3:  low  altitude,  high  dynamic  pressure 

Condition  Ul     medium  altitude,  high  dynamic  pressure 

Condition  -5*  high  altitude,  low  dynamic  pressure 

Condition  6:  high  altitude,  high  dynamic  pressure 
and  are  specifically  referred  to  as  conditions  1  through  6  in  the 
remainder  of  the  report. 

Insertion  of  the  appropriate  transfer  function  into  the  control  loop 
results  in  an  expression  for  the  open  loop  function  H3(s)  which  is,  of 
course,  merely  a  combination  of  the  loop  transfer  functions  into  a  single 
expression 

V8)  °o(s)  °c(a) 


m{s)    =    %  '  e^TJ  aa  '  yrj  Vk  *  G^Tsl  INST. 
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TABLE   I 


Three  Degree  Lateral  Transfer  Functions  for  K_/&u  at  Indicated 

Conditions  of  Flight 

Condition  l; 
Condition  2: 
Condition  3: 

.  a,     -    .g. yi8Ste+o.o<**)(*-4.it>72.X$+s.2iei  tjoAfse) 

Condition  U\ 

N  >c     o.aig^ (s-o.oiZ.eXs+o&LZ£)(s-t.i238Xs+-  I-Q&&) 

Condition  5: 

-       3,3/D4(s-o;fl073X^4.^3Xs--g>%pys-f-4.4gg3) 
V*l      (Sto.0^K6+^77W(5+0.357l±J^.3784.) 

Condition  6; 

H  /&_     0,a*oSY6^.yi)6H-0/?/ffXs-3374lXs+3.4gf) 
r  Tl      (5  +  0,0O7^X5to.'H0+)(6+O./70(,t^/.7OG) 
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which  may  be  examined  by  root  locus  or  frequency  response  techniques  to 
determine  the  linear  characteristics  of  stability  and  response. 

Primary  investigations  of  the  characteristic  responses  were  made  by 
the  root  locus  method.  Since  the  characteristic  equation  of  a  closed  loop 
system  is  given  by 

1  t   m(s)  =  0 
where  the  overall  gain  constant  K  is  the  open  loop  gain  and  G(s)  represents 
the  overall  open  loop  function  in  La  Place  Transform  notation.  It  should 
be  noted  that  the  sign  of  K  may  be  either  positive  or  negative  depending 
upon  the  chosen  positive  response  directions  in  Fig.  1.  The  -  sign 
results  from  the  combination  of  the  sign  of  K  and  the  feedback  sign  used 
to  close  the  loop. 

Open  loop  functions  for  each  flight  condition  are  indicated  in  Table 
II  for  the  case  of  the  aircraft's  lateral  response  being  N Vb-.  The  sub- 
scripts on  each  of  the  functions  refer  to  the  previously  designated  flight 
condition.  In  the  interest  of  brevity  in  all  tables  which  follow  and  list 
lateral  transfer  functions ,  the  denominators  of  the  open  loop  functions 
contained  in  Table  II  shall  be  designated  by  A (s )  where  the  subscript  n 
specifies  the  condition  of  flight. 
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TABLE   II 


Open  Loop  Transfer  functions  for  the  Lateral  Response 
ILAj  at  Indicated  Conditions  of  Flight 


Condition  1; 


m  (B  }1      (±TZ&£fj&^^  7) 


Condition  2: 


1 

Condition  3: 

Condition  k\ 


Condition  5: 


U,5     *(s+o.o44Xs±4.77ilX5+0. 3S7I  t^  Z.  Sie4X<>Ho)(s+V>±j2jo)(sigyt  ftj&j) 


Condition  6: 


KD(8)      ^i.<MO-*iJfa(*-cu>oi)(^-1itsXs*£mX£+S.4Bl) 
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Preliminary  root  locus  plots  were  made  on  these  functions  in  Table 
II  and  they  were  found  to  be  unsuitable  for  both  positive  and  negative 
feedback.  For  positive  feedback  the  locus  emanating  from  the  complex 
poles  immediately  entered  the  right  half  plane ,  or ,  the  unstable  region  of 
the  S-plane.  The  situation  for  negative  feedback  produced  stability  as 
concerned  the  locus  leaving  the  complex  poles.   However,  for  o&  flight 
condition,  a  real  root  was  present,  for  all  values  of  loop  gain,  in  the 
right  half  plane .  The  root  locus  sketches  are  not  contained  in  this  report 
but  mentioned  to  indicate  that  the  use  of  lateral  acceleration  as  the 
principal  dynamic  response  required  compensation . 

At  this  point  it  is  interesting  to  note  the  aircraft  function  poles 
and  zeros  since  these  exist  in  the  critical  portion  of  the  S-plane  and  are 
predominant  in  determining  the  closed  loop  characteristics  of  the  system. 
Fig.  5  indicates  the  characteristic  arrangement  of  the  airframe  poles  and 
zeros  for  a  supersonic  aircraft.  These  differ  from  the  conventional  air- 
craft in  that  the  latter  case  usually  exhibits  two  complex  pole  pairs 0  The 
complex  pole  pair  for  the  supersonic  vehicle  is  indicated  as  the  short 
period  or  dutch  roll  pair  since  the  locus  emerging  from  this  set  usually 
are  dominant  in  fixing  the  short  period,  higher  frequency  oscillation. 
The  real  pair  occasionally  give  rise  to  a  complex  locus  resulting  in  a 
low  frequency,  long  period  oscillation  and  are  referred  to  as  the  phugoid 
or  long  period  pair.  Characteristically,  roots  contained  in  the  phygoid 
locus  are  lightly  damped  while  short  period  roots  may  attain  much  more 
favorable  damping  ratios.  Characteristically,  real  roots  may  also  exist. 
Usually,  one  real  root  is  contained  in  the  negative  real  axle  and  one  real 
root  occurs  about  the  origin  either  in  the  positive  or  negative  real  axis. 
The  former  root  locus  segment  gives  rise  to  a  root  called  the  roll  damping 

factor  which  controls  the  stiffness  of  response  in  roll  to  a  rudder  deflection. 
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The  root  about  the  origin  determines  the  spiral  characteristics  of  the 
aircraft.  As  long  as  the  residue  or  coefficient  of  this  particular  root 
is  small ,  and  should  the  root  exist  in  the  unstable  half -plane  t  the  spiral 
instability  is  easily  compensated  for  by  the  pilot  or  autopilot.  Should 
this  residue  be  significantly  large,  the  aircraft  is  obriously  spirally 
unstable.  Referring  to  the  list  of  transfer  functions  for  N  /bu  in 
Table  I,  it  may  be  noted  that  these  poles  and  zeros  generally  conform  to 
the  above  cases .  The  zeros  may  appear  as  complex  but  generally  as  real 
points . 

Since  compensation  was  desirable  the  suggested  scheme  derived  from 
basic  design  research  conducted  by  Autonetics  design  personnel  was  analyzed. 
This  consisted  of  inserting  a  lead  filter  network  into  the  loop  with  trans- 
fer function 

9~[sl  FILTER  1  ~  (O.ls+1) 

and,  it  may  be  noted,  a  steady  state  gain  of  unity.  By  insertion  of  the 

filter  it  was  desired  that  the  locus  would  be  shaped  so  that  oscillatory 

response  was  determined  largely  by  the  dominant  complex  roots  and  therefore 

second  order  assumptions  could  be  applied  to  regulate  damping  ratio  and 

natural  frequency. 

By  methods  contained  and  illustrated  in  Appendix  A,  solution  of  the 

basic  equation  defining  the  root  locus,  the  characteristic  equation,  (where 

Hj(sL  now  is  defined  as  HG(s)m  •  H3(s)  filter) 
n  n 

i  -  ro(s)n  =  o 

was  carried  out  for  both  conditions  of  feedback.  In  this  instance,  the 
case  of  negative  feedback  proved  to  be  undesirable  since  the  locus  emanating 
from  the  complex  poles  directly  entered  the  unstable  or  right  half  plane. 
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However,  the  condition  of  positive  feedback  produced  loci  which  were  stable 
over  some  variable  gain  region  and  whose  complex  short  period  loci  indicated 
that  a  variable  damping  ratio,  with  varying  loop  gain,  could  be  achieved. 
The  root  loci  plots  of  the  Hy  plus  lead  filter  or  filter  1  response  for 
positive  feedback  are  included  as  Figs.  6  through  11.  Gain  points  on  these 
plots  are  marked  as  values  of  the  variable  loop  gain,  K-  ,  inserted  into 
the  loop  and  not  the  open  loop  gain  as  is  normally  the  case .  Furthermore , 
the  variable  gain,  KL  ,  is  expressed  in  the  units  deg/g  on  all  plots  and 
references  in  this  report.  Since  all  transfer  functions  were  given  in 
terms  of  radian  measure,  KL  was  expressed  as  rad./g  for  all  computations. 
This  convention  was  incorporated  because  of  the  aircraft  industry's  practice 
of  referring  to  control  surface  deflections  in  degree  units. 

In  an  attempt  to  provide  increased  dominance  by  the  complex  roots , 
additional  compensation  was  proposed  by  the  authors  in  the  form  of  an 
additional  lag  pole  in  the  lead  filter  at  £  equal  to  -  0o5o  This  filter 
took  the  form 

°c(s)  .      fa+l) 

0R(s)  FILTER  2     (.ls*l)(2s+l) 

In  addition  it  was  desirous  to  perturbate  the  pole  at  -  10,0  and  the  zero 
at  -  1.0  in  filter  2  to  investigate  the  effects.  The  following  filter  was 
proposed : 

°c(fl)  m  (s+2) 

G~TT5  FILTER  3  ~   (0l25s+l)(2s+l) 

Filters  2  and  3  were  examined  under  conditions  of  negative  and  positive 

feedback  by  the  methods  contained  in  Appendix  10  Root  locus  plots  with 

these  filters  inserted  as  loop  compensators  are  contained  in  Figs .  12 

through  17  for  filter  2  and  Figs.  18  through  23  for  filter  3.  In  each  of 

these  cases  satisfactory  loci  shaping #was  obtained  for  negative  feedback 
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conditions.  Positive  feedback  was  undesirable  since  complex  roots  processed 
directly  into  the  right  half  plane. 

The  final  suggested  implementation  of  N  was  in  combination  with  some 
proportional  amount  of  r,  yaw  acceleration,  resulting  in  a  transfer  function 
(N  ♦  5*  r)/b-.     The  selected  value  of  Ks   determined  the  offset  distance  of 
the  measuring  accelerometer  from  the  e.g.  Positive  K»   indicated  positioning 
of  the  instrument  or  the  center-line  of  the  airframe  forward  of  the  e.g. 
while  negative  K*   placed  the  instrumentation  aft  of  e.g. 

The  open  loop  transfer  functions  for  (N  ♦  ^*  r)/b„  are  listed  in 
Table  III.  According  to  selected  values  of  K>  equal  to  -  0.1,  -  002  and 
0.31  These  values  were  chose  a.  as  indicative  of  the  probable  offset  dist- 
ance of  the  accelerometer  from  the  e.g.  The  variable  loop  gain  was 
designated  as  L  and  expressed  in  the  units  deg./g. 

Root  locus  solutions  were  obtained  using  these  functions  as  open  loop 
transfer  functions.  Due  to  the  enormity  of  the  investigation  of  these 
functions  and  the  fact  that  N  with  compensation  was  being  thoroughly 
analyzed,  studies  ofcly  of  the  feasibility  of  this  conciliation  were  carried 
out.  Family  plots  indicating  the  paths  of  the  loci  in  the  regions  of  the 
origin  and  imaginary  axis  are  contained  in  Figs0  24-  through  30  for  each 
of  the  values  of  Ki  mentioned  above. 
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TABLE   III 


Open  Loop  Lateral  Transfer  Functions  for  (N     ♦  K*  r)/6R 
with  Varying  Values  of  fc  for  Each  Flight  Condition 


K£  s  *o.i; 

Condition  1; 


m,(a)      m    o.\\oz*\o%X^Cs-o-oi&)(s-\.i1Q)(jz+\Joii±ioAitZ^ 


Condition  2: 


Condition  3: 

ro«(s)3  =  ^ ^ # 

Condition  4: 

b,(b)* — "* fiprA 

Condition  5: 

B«(.)5  =  ^ ^ 

Condition  6: 

k        RD»(s)6  =  >)      j^-  
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TABIE   III   (Cont'd) 


*}    =-0.1; 
Condition  1;  1 

m,(s)i  =  * aTSY ' — 

Condition  21 

'2 

Condltltn  3: 

Condition  U\ 

FG'Cs^    = **J ^f- 

Condition  5; 

ffi,(8)5  = ^ zr^T : — 

Condition  6: 

b»(b)6 *«* ^-^ 
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TABIE  III  (Cont'd) 


Condition  1; 


KG 


Condition  2\ 


Ell}2     "  » A(s)l 


Condition  3! 


Condition  4; 

Condition  51 

ro  (sN    =    7      a^_ 


Condition  6: 


K((s)    =  o.4si  x  io*/^  (s-o.o&fX&o.l&sXs-s.aisXs+s.to.) 
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TABIE  III  (Cont'd) 


X±    =    -0.2; 
Condition  1; 

,agaQ 


KG 


.  ,    x  Q.34Q4X  iO*tij„(S'0.6,sXz-0."7SX*+  l-Olb&±<\0:: 

(s)i   =  "• ZTil^ *— 


Condition  2i 


o 


KD«(s)2  -    *** A(a)2 

Condition  3: 

111  (8)3 * *{s)3 

Condition  4: 

-.,  *  _   b.4on  *io*  k>M^fe-Q.OBrTaX^-4.8o35X^^tfoQ(6» 4.S&7) 

k»(b)4  -   =»       a(s)4 

Condition  51 

KD'(s)5  -     "Y         a(s)5 

Condition  6; 

_,,  v  I.(T)3XI0*  l<W  (6-QtooQ^X^-^o^'^fe-2X,y5y5-h2,74?"y) 

m»(8)6  =  =&-> tla)         —  
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TABIE  III  (Cont'd) 


^    =    *o.3; 

Condition  1; 


_  -.0433**0*  «^(s-0.0Sl)(&i>i.*l4X*-i>M(.X 


6+2,/4i) 


Condition  2! 


8B'(8)2     =    «^ ^ 


Condition  3; 


3  A  v3  /o 


Condition  4.; 

E  (b)4  = ^ — airr         — 

4  * 

Condition  5: 

5  °  A(s)_ 

Condition  6: 

A(s)6 
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TABLE  III  (Cont'd) 


K£  =  -0.3 
Condition  1; 

Condition  2; 

ro,(s)     m   l.2tf xuPkUiji [s-aoJi&Ys+2.3t2.X*-i-5*A-X£+i-7S85) 
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3   "  ■ ZT^T3 
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KGi(a)   _  sp ^_l_ 
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In  order  to  corroborate  the  information  obtained  from  the  root  locus 
plots  of  N  and  compensating  filters,  analog  computer  simulation  of  the 
airframe  and  control  loop  was  carried  out  as  detailed  in  Appendix  B. 
The  data  obtained  from  the  analog  study  is  developed  and  discussed  in 
Chapters  on  the  system  analysis,, 
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Chapter  V 

Scope  of  This  Investigation 

Yaw  Rate  ,  r,  as  a  Lateral  Response 

Referring  again  to  Fig.  2,  the  lateral  airframe  response  yaw  rate, 
4.,  may  be  inserted  into  the  control  loop.  The  lateral  transfer  functions 
then  become  r/b~  where  r  and  6p  may  be  expressed  in  deg./sec  or  radians/ 
sec  and  degrees  or  radians  respectively.  The  airframe  equations  for 
each  flight  condition  are  contained  in  Table  IV. 

Combination  of  these  transfer  functions  with  the  instruments  and 
the  variable  gain  parameter,  K  ,  expressed  in  the  units  deg/deg/sec  or 
rad/rad/sec  are  contained  in  Table  V  and  are  the  open  loop  transfer 
functions  where  A  (s )  has  been  previously  defined. 

Root  locus  plots  as  a  function  of  variable  gain,  K  ,  were  obtained 
as  in  previous  examples  and  are  plotted  as  such  in  Figs.  31  through  36. 

As  in  the  preceding  investigations,  it  was  hoped  from  the  anal- 
ysis to  determine  whether  this  lateral  response  offers  the  desirable 
characteristics  upon  which  the  adaptive  scheme  is  to  be  based 0 

Simulation  of  the  aircraft  and  control  loop  were  again  made  on  the 
analog  computer  to  confirm  the  information  derived  from  the  root  locus 
investigation  as  well  as  establish  definite  transient  characteristics. 
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TABIE   17 


Three  Degree  lateral  Transfer  Functions  for  the  Response 
r/bp  for  Each  Specified  Flight  Condition 


Condition  1; 


Condition  2: 


-l.LSOl  (s+o.QOtolt 40,5mgXs+*-254<?) 


r/^  ^ooii^+TimX5-H D.  S31Utj M 172.; 

Condition  3: 

a,  -   -g.3&S66-fl'/3^X^0^7g3;fe*.S\/70 

Condition  4-: 

^  C6-h0.O44)CfH-4.nf2.X^-'-^3S7ltj2,37H; 
Condition  5: 

r/^  (S+O. O  I  44X^0'  3334-XS+  °' '  a7±3f >•  U>T?) 
Condition  6: 

-  a    -2.4-rfgU-0,04,gg  ±4Q.422.l)fe+  I.034&3 

^  (s+OfOOiOte+^+cA  X^-o-no(p±^  met*) 
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TABLE  V 

Open  Loop  Transfer  Functions  for  the  Lateral  Response 

J/ 

'fa   at  Specified  Flight  Conditions 


Condition  1: 


1  A(s)* 

Condition  21 

B.i(8)  3  -2J237XioV/i.^aoofeii^>Sy/6Y5^2,z5^) 


>2     '  ATsT2 

Condition  3:   . 

bii(b)       _  - II . 227€ X/0* *k  fe-ft  /S34 )&Q, 4flf3)fo»5. /7<fc) 

3  A(s)~ 

Condition  4: 

ro , ,  (fl }       m  -ll.gfrl  *lQg&  te^Ptfg*flO/7ag^6+^fo6Q 

4  A(s). 

Condition  5: 

5  a(s)« 
Condition  6: 

ro  (8)6  s iwJT z^ 
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Chapter  VI 
Scope  of  This  Investigation;  The  Adaptive  Scheme 
The  adaptive  scheme  used  to  vary  the  loop  gain  and  achieve  the 
desired  damping  ratio  is  indicated  in  block  diagram  form  in  Fig.  37. 
The  scheme  is  designed  to  select,  according  to  an  error  criterion,  the 
characteristic  damping  ratio  at  the  criterion  minimum  by  furnishing 
through  logic  circuitry  a  voltage  proportional  to  the  error  incurred 
which  in  turn  drives  a  servo  positioned  potentiometer.  Design  of  the 
logic  circuitry  was  accomplished  in  conjunction  with  R.  K.  Smyth,  Senior 
Research  Engineer  and  Supervisor  of  the  Advanced  Analysis  Unit,  Autonetics 
Division,  North  American  Aviation  Corporation,,  # 

The  error  criterion  selected  as  the  basis  of  the  adaptive  scheme 
was 

T 


^  I  A  )  dt 


E  = 

0 

where  A  represents  the  signal  chosen  to  represent  the  particular  airframe 
variable  which  it  is  desired  to  control 0  With  lateral  transfer  functions 
based  on  the  measurement  and  control  of  N  this  variable  would  become  the 
quantity  to  be  minimized.  Furthermore,  since  it  would  be  ideal  as  concerns 
aircraft  control  to  maintain  zero  lateral  acceleration  or  zero  angle  of  yaw 
with  respect  to  the  aircraft* s  course,  a  minimization  of  these  quantities 
would  serve  to  approach  ideal  conditions  with  respect  to  lateral  motion. 


*  Patent  application  filed 
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Reference  5  examines  the  behavior  of  the  integral 

|a|  dt 


$ 


0 

as  a  function  of  damping  ratio  for  a  second  order  servomechanism.  Fig. 
38  illustrates  this  graphically.  While  the  integral  certainly  minimizes 
at  a  damping  ratio  of  0.7  for  the  second  order  servomechanism,  there 
was  no  reason  to  suppose  that  for  a  higher  order  servomechanism,  where 
dominant  complex  roots  occurred,  that  the  criterion  did  not  minimize 
at  the  complex  root  location  resulting  in  the  least  oscillatory  or  the 
most  damped  condition.  Root  locus  studies  detailed  in  preceding  chapters 
indicated  that  at  no  time  would  damping  ratios  greater  than  007  be 
attained „ 

Therefore ,  for  these  reasons  and  the  additional  fact  that  the 
lateral  response  was  to  be  minimized ,  the  error  criterion  stated  above 
was  adopted 0 

The  adaptive  scheme  was  designed  to  operate  by  supplying  periodically 
to  the  airframe  a  disturbance  in  the  form  of  a  short  duration  rudder  pulse . 
The  response  variable  to  be  minimized  was  measured  by  conventional  instru- 
mentation. The  Instrument  output  then  formed  the  basic  signal  operated 
upon  by  the  adaptive  logic  circuits 0 

The  logic  circuitry  is  schematically  illustrated  in  Fig.  39  and 
consists  of  three  memory  storage  elements  or  intergrators  I,  ,  Ia  and  I3, 
relay  switches  and  summation  devices  which  finally  provide  the  output  or 
drive  voltage  to  the  servo  driven  potentiometer 0 

The  time  sequencing  of  the  logic  circuit  was  accomplished  by  relay 
control  of  a  stepping  switch.  The  optimizing  period,  or,  the  period  during 
which  the  logic  circuit  computes  its  output  and  the  control  loop  gain  is 
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adjusted,  is  shown  in  Fig.  40.  The  period  is  divided  into  distinct  intervals 
during  which  the  various  components  perform  the  indicated  tasks.  At  time 
ty,  »  integrators  I,  and  I-  might  integrate  the  value  of  the  criterion  over 
the  period  t-,-1  to  t,-.  I?  does  not  integrate  at  this  time  but  contains 
the  stored  value  of  the  criterion  from  the  previous  interval.  This  occurs 
since  I,  and  I-  alternate,  or,  reset  and  integrate  every  other  period. 
During  period  1,  I,  and  I,  are  the  working  integrators  while  I-  retains  a 
stored  value  from  period  0.  During  period  2,  I?  has  been  reset  to  zero 
and  integrates  with  I,  while  I,  retains  its  integrated  value  from  period  10 

The  criterion  value  as  obtained  by  I-  is  combined  with  some  reference 
value  of  the  criterion  to  form  the  adaptive  error  which  is  fed  to  the 
positioning  servo.  L.  and  I?  compare  values  by  applying  their  voltages  to 
a  relay  switch.  If  I,  >  I«  J  ,  during  ty,  to  t,_  mentioned  above,  the 
sign  switch  reverses  since  the  value  of  error  is  increasing  and  the  sign  of 
the  adaptive  error  applied  to  the  pot  servo  is  reversed .  Thus  the  gain 
will  be  increased  by  an  increment  proportional  to  the  adaptive  error  in  the 
opposite  direction.  Should  L.  \  K.     I?  I  ,  the  sign  will  not  reverse  and  the 
gain  will  vary  in  the  same  direction  as  in  the  previous  period.  Fig.  41 
indicates  the  operation  of  these  integrators  within  several  typical 
optimizing  periods. 

In  order  to  insure  good  response  from  the  servo  driven  potentiometer, 
the  servo  gain  was  made  a  function  of  the  value  of  the  adaptive  error. 
A  threshold  value  was  incorporated  below  which  the  servo  gain  was  a  fraction 
of  its  gain  with  an  above  threshold  value  of  adaptive  error.  This  was 
intended  to  reduce  the  sensitivity  of  the  gain  adjuster  when  the  minimum 
value  of  the  criterion  was  approached  and  thus  reduce  hunting  of  the 
mechanism.  A  non-linear  potentiometer  was  included  also  to  permit  more 
sensitivity  in  gain  selections.  The  pot  characteristics  are  presented 
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la  Fig.  4-2  and  as  indicated  provided  a  lower  rate  of  gain  change  in  the 
lower  gain  regions  as  a  function  of  angular  positioning. 

The  criterion  reference  signal  was  available  as  a  predetermined  DC 
signal  of  desired  polarity  or  was  capable  of  being  generated  as  some 
proportional  increment  of  the  measured  criterion  value  during  the  optimizing 
period.  The  latter  method  was  desirable,  in  case  it  was  necessary,  to  off- 
set noise  effects  in  the  system  by  error  derivations  from  signals  which 
included  the  same  noise  effects  thereby  tending  to  cancel  undesirable 
inputs . 

The  selectidn  of  the  time  interval  durations  within  the  optimizing 
period  was  made  on  the  basis  of_  the  stability  analysis  conducted  on  the 
control  loop  parameters  0  The  primary  requirements  were  that  the  intervals 
allow  integration  over  the  maximum  rate  of  error  change  and  that  sufficient 
time  be  allowed  for  comparison  and  reset  of  integrators  in  the  desired 
sequenof 9     It  was  determined  from  root  locus  plots  of  desirable  lateral 
responses  that  a  choice  of  2  second  time  interval  durations  would  be 
adequate  since  in  this  period  of  time  all  responses  had  achieved  one  over- 
shoot  and  essentially  one  undershoot.  Thus  the  most  sizeable  variations 
in  response  were  measured. 

The  Analysis  and  basic  testing  of  the  adaptive  scheme  was  carried  out 

« 

by  analog  computer  simulation  of  the  three  degree  of  freedom  aircraft  and 
integrating  units  of  the  adaptive  controller.  Timer  functions  were  per- 
formed by  a  mechanical  stepping  switch  which  actuated  portions  of  the 
adaptive  controller  according  to  the  previous  discussion  in  thie  chapter. 
In  order  to  provide  a  certain  amount  of  testing  flexibility,  provision 
was  made  at  the  computer  console  to  provide  a  means  of  altering  the 
magnitude  of  the  rudder  pulse  disturbance  to  the  control  loop  and  to 
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permit  an  unrealistic  instantaneous  switch  from  one  flight  condition  to 

i 

another 0  While  unrealistic,  this  enabled  a  maximum  test  of  the  adaptive 
controllers  ability. 

A  limited  investigation  of  the  adaptive  scheme  was  performed  on  the 
analog  computer  by  applying  the  controller  to  a  simulated  second  order 
servomechanism  whose  maximum  damping  ratio  was  adjustable.  The  sole 
purpose  of  this  investigation  was  to  permit  some  check  on  the  character- 
istics of  the  adaptive  controller  as  applied  to  a  basic  system  prior  to 
its  application  to  the  more  complex  aircraft.  It  was  determined  that 
the  adaptive  logic  was  correct  and  that  whenever  the  maximum  damping  ratio 
was  0.7  or  less,  the  logic  was  eapable  of  adjusting  the  servomechanism 
gain  so  that  this  maximum  damping  ratio  was  achieved. 


76 


CHAPTER  VII 
The  Analysis  of  Ny  Plus  Lead  Filter  System 

As  stated  previously  in  Chapter  IV,  the  case  of  negative  feedback 
using  filter  1  was  entirely  unsatisfactory.  Using  negative  feedback  a 
damping  ratio  greater  than  that  for  IC.  equal  to  zero  could  not  be 
obtained.  Therefore,  negative  feedback  will  be  discussed  no  further  for 
this  case . 

The  case  of  positive  feedback,  however,  shows  satisfactory  results 
possible  in  several  ways.  First,  it  permits  adaptive  control  of  the 
feedback  pain.  By  this  it  is  meant  that  an  improvement  in  damping  ratio 
can  be  obtained  by  varying  IC.  „  Secondly,  the  period  of  the  high  frequency 
pair  is  acceptable,  and  finally,  the  range  of  values  for  K„  is  reasonable. 
Table  VI  contains  the  estimated  optimum  values  of  K„»  damping  ratio,  and 
natural  frequency,  oo  ,  for  this  scheme.  To  obtain  these  values  the  system 
was  considered  to  approximate  a  second  order  system.  The  value  of  K- 
was  chosen  from  the  root  locus  plots  such  as  to  give  the  best  damping  ratio 
and  still  have  stable  roots  of  the  characteristic  equation.  From  these 
values  one  can  see  that  a  mean  £  of  0.3  is  achievable,  excepting  condition 
3.  In  other  flight  conditions,  2  and  6,  greater  damping  ratios  can  be 
obtained.  The  period  of  the  short  period  response  varies  from  3.5  seconds 
to  1.36  seconds. 

In  figure  6  through  11,  and  in  all  subsequent  root  locus  plots,  no 
attempt  was  made  to  show  the  loci  emanating  from  the  autopilot  actuator 
and  measuring  instrument  poles.  This  was  possible  since  the  IBM  studies 
showed  that  these  roots  moved  only  a  short  distance  for  the  range  of  loop 
gain  employed  and  the  roots  could  be  considered  as  essentially  constant. 

Table  VII  contains  a  tabulation  of  several  of  the  significant  root 
loci  characteristics  for  all  six  flight  conditions  of  this  case.  Note 
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TABLE  VI 

Root  Locus  Second  Order  Approximations  of  Transient 

Characteristics 


N     Plus  Filter  1 
7 


Flight 
Cond. 

deg/g 
20 
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rad/sec 
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0o312 

1.8 
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0M 

202 

3 

0.6 

0o22 

406 

k 

0.7 

0.29 

3o8     . 

5 

20 

0.25 

2.0 

6 
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0.42 

20l 

78 


1 

9 


I 

to 
o 

«H 
P 

to 

•H 
P 

$ 
V 

Crf 
p 

crf 

43 
o 

(0 
O 


P 

o 

o 

CD 


o 
o 


Pi 

■H 

•H 


a 

p 

3 

s 

a) 

<H 

ia 

CD 

s 

CO 

o 

•H 

P 

CO 

«H 

p 

Pi 

o 

CD 

o 

P 

cr; 

o 

erf 

r-\ 

p 

crf 

a 

£ 

o 

CO 

o 

•H 

•8 

P 

CO 

•H 

•H 

P 

Pi 

a) 

CD 

D- 

P 

O 

p 

Crf 

P 

Pi 

o 

Crf 

43 

X 

00 

O 

ID 

o 

•H 

p 

CO 

•H 

P 

0) 

73 

p 

•H 

o 

O 

IP 

2 

erf 

£ 

43 
o 

• 

• 

«M 

bO 

& 

•H 

P 

• 

■a 

•a 

•H 

o 

rH 

o 

t^ 

T3 

CO 

•H 

■g 

CD 
P 

X 
erf 

c 

a 

■H 

CD 

•H 

erf 

p 

P 

rH 

rH 

8. 

C 

erf 

T3 

& 

0 

£ 

s 

CD 

■e  s 

• 

•H 

rH 

43 

o 

o 

bO 

$ 

X 

o 

8 

<D 

CO 

H 

c 

CO 

CD 

rH 

■ 

0 
O 

a 

a 

P 

T3 

Cm 

CD 

• 

•H 

43 

to 
• 

8 

3 
O 

EC  rH  TJ 

CD    O    P, 

TJ  «m    CD 

rH 

•«H  43 

t*«  o  3 

•  o  o 

P3  rH  73 


•a 

CD 

S«r>  «H 
O  Cm 
Crf  «H   «H  «H 

erf  p  -a 

fe  bO   CD     O 
rH  P,  6 

rQ    P 
•H    CD    P    W 

co  x:  p  p 

10  wo   o 

0-H£    O 

Pi  43    co   p 


P 

crf 

CO 

3 
o 
o 


I    I 

Pi     C 

It 

P   o 

O  TJ 


S. 


rH       o    S 

O    P. 

• 

•      o   o 

o 

tt,     O    «M 

CO 

•rH 

o 

CC     1     CD 

CO 

p 

fell    -H 

£     CO     U 

d 

crf 

>£> 


CO  Cm  rH 


a 


CD    O     CX 

3    P 

TJ  Cm    q> 

rH 

•■H  43 

to  o3 
•  o  o 


p 

Pi 


C  rH 

•h  a 


p 

a 


o 

T5 

I 
P 
O 

to 


Cm 
H 

crf 


Pi 

■H    0) 

as 

•H 

•    CO 

O     CO 

CO     O 

O  Pi 


a 

3  S 

<m  • 

■H  En 


8 

O 
55 


CV 


P>    Pi 
P    O  «H 


bo  a, 
p 

CD 
> 


s 

o 
p 

Cm 


P 


TJ  >-i 


P 

a 

p 
p 

o 

43 

CO 


J 


CD     P< 


■P 
P 

O    CD 
b0>-3    P.H 

CD    Crf 

P    CE  rd  P 

CO 


-j  C  43 

P  crf    3 

C  rH     O    C 

CD  P.TJ  *H 


£3  O 
CD  O 
bO  P 
P 

CD  l~t 
>    crf 


2 


fe 
S 


S  «H 
feg 


co 


P 

8 

p 

CD 
<-i 
XI 

Crf 

p 

CO 


O  • 

P      •  CD 

C    CD     I    rH  O 

•H    C    CO    Crf  C 

»0           Crf  «H    CD  CD 

O    CD  rH  TJ    P  bO 

•H    >•    P,  P 

*"    fa  .       ®    fe»  P 

CD    7  Cm    Q  ju  > 

P.  o  H    a  «H 

erf    erf  TJ  TJ 
■P    CO  X    C    CD 

U  P          «H  P  +> 

o  o    •  e  P  Q 

J3    O  P    O    O  5 

CO    P    P  Q  P  P 


8 


CO 


c> 


p  p 

$% 

P    CD     I 

erf   >    C 
P.  O  «H 

rH         P 
crf     *  O 

P  <^  o 

CD  P 

C      • 
CD  'Q  rH 
bO  C    Crf 

g8£ 

Sao 

CO    crf    c 


TJ 
O 


t 

P 
P    P 

O  «H 
X  Crf 
CO    P. 


§CD  Xj 
CD  rH    O 
O  rH  43    erf 

2  H+,  i 

CD  t-i 

•rH     C  > 

>»43  -H  O       • 

■P    3  rH     O 

■HOB  CO 

H-O    O  *  O 

■H           «H  CD 

43    S    crf  C  43 

erf    O    E  erf    bO 

P    P    CD  r-\  «H 

CO  Cm    P  P.43 


In  rt 
plane 

CD  «m 
C  rH 
O    crf 
Z  43 

a 


P  P  «H 

P     P     Crf     Crf 


£ 


O 


79 


1 

o 


9 


3 

ct? 
-P    <D 

&  3. 


4> 

n 

•H 

J-P    »h 
O 

r-i     U 

ce  3 

is 


01 

o 

O     01 
•H 
fn     M 

V 

■P   ocf 

o  3 

X!  X! 

CO  o 


«9 

O 
•H 
■P 

CO 
«H 

h 

•O  -P 
•H  O 
O    o5 

3  05 
»P  •*= 


Cm  t)C 
CO  -H 
«   fc 

■S.'S 

•H    O 


a 

i 

§ 

patter 
sted  as 

• 

fe 

12    *H 

»H 

o 

I   M 

<H 

3  a) 

• 

CO 

CO 
•H 

-a 

o5 

13 

o 

1 

£ 

•  <H 

o 

rH    O 

Cm 

«^ 

s 

§a 

fa 

•o  <► 

• 

•h2 

• 

8§ 

WH-O 

£5 


TJ    CO 

S*TJ    0)  +* 
O  «H    O 


etf  «H  «H  Cm 


» 


05  k  «h  t-,  xj  s:  «H 

r-\  5,  O      •   O    O    O 

<£•  fa        B  fc,  tJ       T> 

•H    <D  4>  •  ♦  • 

OlXifafttCGcpcPfa 
01  15)  O  «H  OiHCO-3 

O  «H  X!    aj    ►»  fa    Q«H    gj 


£ 


*s 

0 

0 

fa 

O 

Cm 

CO 

a5 

i 

05 

• 

CO 

3 

3 

8 

0 

rH 

T3 

• 

B 

(I) 

fr« 

O  rH 

• 

b 

O 

ttl    «M 

a 

M 

s 

05 

a  r- 

«H 

p. 

8 

Cm 
rH 

0 

05 

m 

X! 

»T\ 


& 

«M 

01  v> 

05 

i'a 

3  o 
co  o 


8 

o 


vO 


80 


that  no  phugoid  oscillations  will  occur,  and  only  condition  3  has  a  spiral 
divergency  condition.  The  system  exhibits  the  possibility  in  conditions 
1,  2,  5t  and  6  of  a  higher  frequency  mode  going  unstable  first  if  K-  is 
set  at  too  high  a  value.  These  are  also  the  conditions  which  are  relative- 
ly insensitive  to  gain  changes,  especially  with  regard  to  the  movement  of 
the  short  period  pair  of  roots .  Flight  conditions  3  and  U   represent  the 
most  sensitive  conditions  of  the  system.  Here  the  value  of  IL,  must  be 
low  and  closely  regulated  for  stability.  The  varying  degrees  of  sensi- 
tivity to  gain  changes  are  due  to  the  airframe  transfer  functions  and  must 
be  tolerated. 

It  may  be  noted  that  from  these  root  loci  studies  an  acceptable 
basis  for  the  proposed  adaptive  system  may  be  found.  A  variation  of  JL. 
as  flight  conditions  vary  will  produce  a  system  which  will  have  a  more 
acceptable  response  than  a  fixed  gain  system.  This  fact  is  substantiated 
by  the  variance  in  JL.    necessary  to  produce  acceptable  response  over  the 
region  of  flight.  If  the  relative  degrees  of  sensitivity  to  gain  changes 
does  not  prove  to  be  a  prohibitive  factor,  this  scheme  is  suitable  for 
adaptive  control  using  only  gain  changes , 
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CHAPTER  VIII 
The  Analysis  of  Ny  Plus  Lag-Lead  Filter  System 

Negative  feedback  was  used  in  this  scheme  with  the  most  success. 
With  positive  feedback  the  short  period  pair  went  directly  into  the 
right  half  plane,  and  therefore  will  be  discussed  no  further  for  this 
case. 

The  scheme  of  using  negative  feedback  with  filter  2  will  be  dis- 
cussed first.  This  scheme  also  showed  desirable  characteristics  since 
it  was  possible  to  vary  the  location  of  the  dominant  complex  pair  by  a 
variation  of  YL.   .  However  the  frequency  associated  with  this  pair  of 
roots  was  in  general  lower  than  those  experienced  with  filter  1  and  the 
gain  variation  range  was  considerably  larger.  At  the  same  time,  the 
damping  ratios  were  generally  greater  and  showed  less  variation  from  one 
flight  condition  to  another  at  the  optimum  values  of  K-  .  The  most 
desirable  values  of  Kw_»Cf  »  and  ca  are  tabulated  in  Table  VIII.  These 
values  were  obtained  by  overall  stability  considerations  and  the  use  of 
a  second  order  approximation  as  in  the  case  of  filter  1.  The  value  of 
IC.  was  again  chosen  for  the  best  damping  ratio  consistant  with  stability. 
The  values  for  this  scheme  using  filter  3  are  not  included  in  Table  VIII0 

Table  VIII  shows  that  a  damping  ratio  of  0o4.  or  Jbetter  is  a  capability 
at  each  flight  condition,  with  0.54  being  the  maximum.  The  period  of  the 
high  frequency  pair  varies  from  4.8  seconds  to  2.1  seconds.  However,  in 
this  scheme,  the  period  is  more  consistent  than  with  filter  1.   In  addition 
there  is  a  capability  of  providing  a  constant  damping  ratio. 

Several  of  the  significant  characteristics  of  the  root  loci  for  the 
six  flight  conditions  have  been  tabulated  in  Table  IX.  This  tabulation 
is  the  case  of  filter  2  only.  The  pugoid  pair  of  roots,  if  oscillatory, 
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TABLE  VIII 

Root  Locus  Second  Order  Approximations  of  Transient 

Characteristics 
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have  a  very  low  frequency  and  move  rapidly  into  the  right  half  plana  once 
they  break  out  of  the  negative  real  axis .  This  results  in  these  roots 
generally  dictating  the  limiting  values  of  IL,  for  stability.  The  pair 
of  roots  on  the  loci  from  the  double  pole  at  a  equals  -10.0  do  threaten, 
in  several  flight  conditions,  to  go  unstable  at  about  the  same  values  of 

The  root  loci  study-  using  filter  3  to  compensate  the  lateral  charact- 
eristic Ny  was  conducted  primarily  as  an  investigation  of  filter  perturb- 
ations .  Again,  as  with  filter  2,  positive  feedback  was  unsatisfactory. 
It  failed  to  allow  an  increase  of  damping  ratio  for  an  increase  in  K„  due 
to  the  fact  that  the  short  period  loci  moved  immediately  into  the  unstable 
right  half  plane.  Referring  to  Figs.  18  through  23,  note  that  the  loci 
followed  very  closely  the  characteristic  loci  for  filter  2  at  each 
condition  of  flight.   There  are,  however,  two  major  differences.  First, 
the  loci  from  the  two  real  poles  (at  s  equal  to  -8.0  and  -10 „0)  emerges 
at  a  value  of  s  between  -P .9  and  -9.0.  Secondly,  the  roots  move  along 
the  locus  segments  very  slowly  and  consequently  a  change  in  IC.  produces 
less  movement  with  filter  3  than  it  did  with  filter  2.  These  loci  show 
that  an  adaptive  control  might  be  feasible  but  this  filter  was  not  con- 
sidered further  in  the  study  because  of  the  following  reasons:   (1) 
Higher  values  of  gain  could  be  used  to  get  an  acceptable  damping  ratio 
but  this  in  itself  is  not  an  advantage  over  filter  2.   (2)  The  higher 
gains  which  are  necessary  to  obtain  comparable  values  of  damping  ra,tio 
and  natural  frequency  cause  the  phugoid  pair  to  become  unstable . 

Therefore,  the  scheme  using  filter  3  showed  no  great  improvement 
over  the  advantages  of  filter  2. 

The  study  of  filter  2  and  3  indicates  that  system  characteristics 
desirous  for  adaptive  control  are  present.  As  previously  stated,  filter 
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2  is  preferable  to  filter  3.  While  greater  damping  ratios  ire  obtainable 
with  filter  2  than  those  obtained  in  the  scheme  using  filter  1  the 
corresponding  frequency  is  lower.  The  range  of  YL.     necessary  for  adaptive 
control  using  this  scheme  is  large  and  might  be  difficult  to  obtain  satis- 
factorily by  adaptive  control.  This  range  extends  from  a  ?C.  of  1.4.  to 
130  deg/g  to  insure  provision  for  the  most  acceptable  response.  Great 
accuracy  and  sensitivity  in  the  change  of  fain  by  the  adaptive  controller 
would  be  required  since  the  stability  limit  of  several  of  the  flight 
conditions  is  close  to  the  most  acceptable  £ain  ror  short  period  root 
location.  The  gain  change,  from  one  flight  condition  to  another,  must 
be  quickly  achieved  in  order  to  prevent  large  gain  change  requirements 
from  causing  temporary  instability. 
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CHAPTER  IX 
The  Analysis  of  Ny  Plus  r  System 

The  root  locus  study  of  this  scheme  was  performed  in  the  manner 
indicated  in  Appendix  A.  Although  both  positive  and  negative  feedback 
were  developed,  the  use  of  positive  feedback  was  found  to  be  unsatisfactory 
in  comparison  with  the  loci  using  negative  feedback.  This  again  was  due  to 
the  short  period  pair  moving  directly  into  the  right  half  plane .  Conse- 
quently damping  ratios  could  not  be  obtained  which  were  greater  than 
those  achieved  with  K„  equal  to  zero.  Therefore,  the  case  of  positive 
feedback  is  not  further  developed 0 

Since  the  nature  of  this  ^particular  study  was  merely  to  establish 
feasibility,  only  the  loci  of  the  short  period  pair  was  plotted.  However, 
the  phugoid  pair  loci  was  plotted  for  one  flight  condition  as  a  general 
indication  of  performance.  Compensation  was  not  used,  and  in  general  the 
root  loci  allow  a  satisfactory  placement  of  the  short  period  pair  of  roots 
without  compensation.  These  root  loci  plots  are  a  family  of  loci  in  the 
sense  that  they  show  the  effect  of  a  variation  of  FO  .  For  the  other 
flight  conditions  the  phugoid  pair  either  did  not  become  complex  or,  for 
the  range  of  gains  used,  did  not  become  unstable. 

By  considering  that  the  system  could  b»  approximated  as  a  second 
order  system,  and  choosing  a  value  of  Kl  equal  to  ♦  002  as  representative, 
the  values  of  damping  ratio,  natural  frequency,  and  FL,  used  to  position  the 
roots  are  tabulated  in  Table  X.  These  values  are  again  chosen  so  as  to 
give  a  good  damping  ratio  and  still  have  a  stable  system.  The  range  of 
fC,  is  about  115  to  1  which  might  prove  difficult  to  provide  by  adaptive 
means.  Another  complication  is  that  the  flight  conditions  which  require 
the  smaller  values  of  IL.  are  the  most  sensitive  and  cannot  tolerate  large 
gains  without  severe  instability.  The  available  frequency  associated  with 
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TABLE  I 

Root  Locus  Second  Order  Approximations  of  Transient 

Character is  tics 


N  with  K£  r  ;  K^  =  ♦  0o2 


Flight 
Cond. 

deg/ar 
230 

0.63 

rad/sec 

1 

1042 

2 

15 

0.465 

0.96 

3 

502 

0.575 

2018 

U 

1.7 

0.62 

1.06 

5 

230 

0.476 

1026 

6 

11.5 

0.656 

lo08 
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the  short  period  roots  is  low  and  not  too  consistent.  This  short  period 
varies  over  a  range  of  6.5  to  209  seconds. 

Perhaps  a  brief  look  at  the  root  locus  plots  of  the  individual  flight 
conditions  would  be  in  order.  In  Fig.  24 »  for  flight  condition  1,  the 
short  period  pair  displays  a  varying  behavior.  With  larger  negative  values 
of  K»   the  loci  parallels  the  imaginary  axis  toward  the  origin.  However 
with  smaller  values  of  positive  Kz   the  loci  approximately  parallels  the 
negative  real  axis.  Finally  with  K^  equal  to  0o3  the  loci  starts  to  curve 
up  quickly.  The  effect  of  varying  K>  is  quite  pronounced  in  this  flight 
condition. 

For  flight  condition  2,  as  shown  in  Fig.  25,  the  effect  of  K»  is  much 
less  pronounced.  While  the  loci  does  show  some  shift  due  to  varying  the 
value  of  FOj,  in  general  they  proceed  downward  and  nearly  parallel  to  the 
imaginary  axis.  Ultimately  they  proceed  down  to  the  negative  real  axis  and 
become  real  roots  0 

Fig.  26,  that  of  condition  3,  the  same  pattern  is  shown  as  for 
condition  2.  However,  the  loci  bow  out  more  from  the  imaginary  axis  and 
the  effect  of  Kl  is  more  pronounced.  This  flight  condition,  as  always, 
exhibits  sensitiveness  of  FC,  changes.  None  of  the  loci  curve  away,  rather 
all  of  them  proceed  down  to  the  negative  real  axis . 

Flight  condition  4,  shown  in  Fig.  27,  follows  closely  the  general 
pattern  found  in  the  condition  2  loci.  The  flight  condition  is  sensitive 
to  gain  changes  as  was  to  be  expected.  The  effect  of  K*  was  noticeable  but 
again  produced  no  marked  effect. 

Flight  condition  5,  normally  a  well  behaved  condition,  was  again  so0 

This  is  noted  in  Fig.  28,  where  the  loci  return  to  the  same  general  pattern 

as  that  found  for  condition  1.  The  effect  of  K°  is  shown  by  the  tendency 

of  negative  values  to  pull  the  loci  close  to  the  imaginary  axis ,  Positive 
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values  tended  to  move  the  loci  out  into  the  left  half  plane  and  more 
closely  parallel  with  the  negative  real  axis.  Fig.  29  is  a  plot  of 
the  phugoid  pair  for  flight  condition  5»  to  an  expanded  scale.  This  is 
included  to  show  the  general  pattern  of  this  pair  of  roots  in  all  flight 
conditions  as  well  as  its  pattern  for  condition  5.  The  effect  of  10  is 
not  very  large.  The  roots  follow  pretty  much  the  same  path,  or  at  least 
form  a  very  narrow  pencil.  The  phugoid  pair,  if  K„  is  large  enough  to 
permit  it  to  develop,  contribute  a  very  low  frequency.   It  is  also  seen 
that  the  value  of  gain  chosen  in  Table  X  must  be  reduced  with  positive 
values  of  K£  to  have  a  stable  phugoid  pair.  This  reduction  would  be  less 
than  10  io    of  the  tabulated  val'ue  0 

Finally,  flight  condition  6  is  portrayed  in  Fig.  30.  The  loci  pattern 
here  is  somewhere  between  that  of  conditions  2,  3,  and  U   and  that  of  con- 
ditions 1  and  5.  The  effect  of  Ki   is  noticeable  but  not  pronounced.  This 
loci  family  is  not  very  sensitive  to  gain  changes  and  for  a  range  of  IL. 
values  from  10  to  15  exhibits  almost  constant  frequency  and  only  varies 
damping  ratio.  The  variation  of  damping  ratio  is  quite  small. 

In  summary  it  may  be  noted  that  from  the  standpoint  of  placement  of 
the  short  period  pair  and  stability,  this  scheme  could  provide  some  basis 
for  self  adaptive  control.  It,  of  the  schemes  previously  analyzed, 
provides  for  the  most  consistent  damping  ratio  which  is  closest  to  0.7. 
However,  achievement  of  this  damping  ratio  for  all  flight  conditions 
requires  a  wider  variation  of  JL.     than  the  scheme  analyzed  in  Chapter  VIl£ 
Also  the  frequency  of  the  system  is  low  and  speed  of  response  would  con- 
sequently suffer. 

Perhaps  the  most  interesting  thing  about  this  scheme  is  the  effect 
of  K>  .  The  value  of  K°  has  practically  no  effect  on  damping  ratio  for  a 
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fixed  value  of  1L.   .  Of  course,  one  must  restrict  this  statement  to  the 

%t  ny 

more  well  behaved  regions  of  the  loci0  The  sole  effect  of  K£  seems  to 
be  to  vary  the  frequency  of  the  short  period  pair.  The  dominant  frequency 
increases  with  increasing  values  of  K»   While  it  affects  the  other  roots 
its  effect  is  most  marked  on  the  short  period  pair.  This  particular 
characteristic  may  be  noted  in  flight  condition  5»  Fig.  28.  Here  for  a 
fixed  value  of  KL.  f  the  short  period  roots  for  various  value*  of  K£  lie 
very  close  to  a  constant  damping  locus 0 

There  are  two  more  noticeable  characteristics  of  this  scheme.  First, 
the  lines  of  constant  K„  have  a  negative  slope  for  1L.    values  low  in  the 
stability  range 0  As  IL.  increases  toward  the  maximum  stable  gain,  this 
slope  of  the  constant  K„  lines  becomes  positive  0  This  can  be  most  readily 
seen  in  Fig.  26.  This  shows  that  for  higher  values  of  K>  the  movement  of 
the  root  for  a  change  in  K„  is  increased 0  Second,  the  effect  of  varying 
10  has  the  opposite  effect  on  the  phugoid  pair  as  compared  to  the  short 
period  pair.   It  may  be  noted  in  Fig,  29  that,  as  K>  increases,  the  loci 
tends  to  collapse  on  the  origin.  Thus,  the  frequency  of  this  pair  of  roots 
decreases  with  increasing  K*0 
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CHAPTER  X 
The  Analysis  of  the  Yaw  Rate  System 

The  origin  of  the  root  loci  plots  for  this  scheme  were  obtained 
using  the  procedures  of  Appendix  A.  While  the  loci  for  both  negative 
and  positive  feedback  were  computed,  the  negative  case  was  unsatis- 
factory. For  this  case  the  short  period  pair  did  not  possess  the  desired 
characteristics  and  there  was  a  real  root  at  all  times  on  the  positive 
real  axis  whJch  moved  rapidly  away  from  the  origin  with  increasing  gain. 
Further  discussion  of  this  feedback  case  is  considered  unnecessary. 

The  loci  of  the  control  loop  with  positive  feedback  are  plotted  in 
Fig.  31  through  36  for  each  of  the  six  flight  conditions.  The  values  of 
damping  ratio,  natural  frequency,  and  variable  loop  gain,  K  ,  are  listed 
in  Table  XI.  This  again  assumes  that  the  transient  oscillation  charac- 

* 

teristics  can  be  approximated  by  the  second  order  short  period  pair. 
These  values  were  selected  for  the  damping  ratio  closest  to  0o7  with  the 
loop  still  stable. 

Noting  the  values  in  Table  XI  it  can  be  seen  that  this  scheme,  as 
far  as  the  root  loci  is  concerned,  is  capable  of  providing  desirable 
oscillating  characteristics  throughout  the  region  of  flight  by  a  variation 
in  K  .  The  range  of  gain  is  of  an  order  of  magnitude  readily  attainable 0 
A  damping  ratio  of  0.7  is  possible  with  the  exception  of  conditions  3  and 
J+,     For  these  conditions,  the  most  acceptable  damping  ratio  is  lower. 
For  most  flight  conditions  the  frequency  is  low  but  practically  uniform, 
while  the  frequency  for  conditions  3  and  U   is  high.  An  oi  of  2.0  rad/sec 
is  slow  and  one  of  A, 4.  rad/sec,  as  in  condition  3»  is  slightly  high.  It 
was  noted  that  the  flight  condition  with  higher  short  period  frequencies 
also  exhibited  smaller  damping  ratios.  This  results  in  a  marked  difference 
in  responses  over  the  flight  region. 
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TABLE  XI 
Root  Locus   Second  Order   Approximations  of  Transient 

Characteristics 

Yaw  Rate 


Flight 
Cond. 

Kr 
deR/k 

1.8 

0o7 

rad/sec 

1 

2o0 

2 

0.9 

0o7 

109 

3 

0,2 

0o35 

kj> 

L 

0.3 

0.61 

4.1 

5 

2o0 

0.7 

2025 

6 

0o75 

0.7 

2.U 
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Table  XII  contains  the  more  significant  characteristics  of  the  root 
loci  for  this  case.  Here  again  the  stability  limit  for  Kj.  is  set  by  the 
phugoid  pair  o^roots,  except  for  condition  3.  This  condition  has  a 
spiral  divergency  condition  at  higher  gains  and  this  root  thus  limits  the 
value  of  KJ.  which  can  be  accepted. 

In  summary,  this  scheme  indicates  that  a  variation  in  control  loop 
gain  could  provide  desirable  damping  ratios  and  frequencies  from  the 
short  period  pair.  The  range  of  K  is  small  and  capable  of  provision  by 
a  self  adjusting  system.  However,  the  frequency  of  the  short  period  pair 
is  too  low  in  general  based  on  human  pilot  preferences  previously  indi- 
cated and  varies  by  a  factor  of  2  within  the  six  flight  conditions. 

The  fixed  gain  performance  of  this  scheme  is  quite  acceptable. 
However,  the  response  could  be  improved  with  self  adjusting  control. 

The  phugoid  pair  was  not  troublesome,  but  as  previously  stated,  one 

flight  condition  did  exhibit  spiral  divergence.  Further  discussion  of 

r 

this  scheme  and  its  acceptability  is  contained  in  Chapter  XI, 
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CHAPTER  XI 
Summary  of  Analysis  and  Analog  Computer  Verification 

As  mentioned  in  the  previous  chapters  the  analysis  of  the  system  was 
specifically  directed  toward  noting  the  characteristics  of  damping  ratio 
and  frequency  associated  with  the  dominant  complex  roots.  The  existence 
of  a  phugoid  pair  of  roots  and,  or,  real  roots  close  to  the  origin  will 
certainly  modify  the  transient  response.  However,  in  that  phase  of  the 
analysis  it  was  feasible  to  consider  only  the  phugoid  pair  and  any  real 
roots  from  a  qualitative  standpoint.  This  allowed  the  various  schemes 
to  be  compared  from  the  standpoint  of  root  locations,  since  the  primary 
desire  was  to  place  the  short  period  pair  at  a  desirable  location  for 
oscillatory  dominance . 

The  root  loci  in  the  analysis  took  into  account,  computation,  all 
poles  and  zeros.  However,  only  one  quadrant  and  the  loci  most  important 
to  the  analysis  were  shown  in  the  figures.  Thus  the  analysis  was  based 
on  very  accurate  root  loci  configurations  which  enabled  good  qualitative 
comparison . 

In  this  analysis ,  which  of  course  was  to  outline  the  requirements  and 
feasibility  of  self  adapting  control  by  gain  variation  only,  the  values 
of  gain  considered  best  for  each  flight  condition  have  been  chosen  and 
tabulated  in  Tables  VI,  VIII,  X,  and  XI.  These  were  determined  on  the 
basis  of  a  damping  ratio,  closest  to  0.7,  with  the  system  still  stable, 
and  a  frequency  as  close  to  3  rad/sec  as  possible.  These  were  considered 
to  be  the  ideal  response  characteristics  for  the  dominant  short  period 
pair.  A  reasonable  contribution  from  a  phugoid  pair  could  be  tolerated, 
as  could  a  spirally  divergent  condition,  as  long  as  they  were  within  the 
control  capabilities  of  the  pilot  or  autopilot. 
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From  basic  research  and  previous  experience,  there  exists  the 
definite  possibility  that  the  small  value  real  roots,  very  close  to  the 
origin,  will  in  effect  vanigh  when  the  entire  autopilot  and  all' its  loops 
are  closed. 

As  previously  stated  the  purpose  of  this  analysis  was  to  choose  a 
scheme  which  would  be  the  best  for  self  adaptive  control.  The  scheme 
chosen  was  N  with  filter  1  compensation.  This  scheme  was  chosen  since 
it  offered  the  best  combinations  of  aerodynamic  suitability  and  root 
location  possibilities. 

As  this  rudder  control  loop  normally  functions  independent  of  the 
maneuvering  section  of  the  autopilot,  this  loop  plays  a  dominant  role 
in  insuring  a  coordinated  turn  for  the  aircraft.  This  coordination  is 
mandatory  for  acceptable  maneuvering  performance .  Since  lateral  accel- 
eration of  the  e.g.  is  one  of  the  best  indications  of  coordinated  turning 
there  are  no  aerodynamic  reasons  for  rejection.  With  N  equal  to  zero  the 
longitudinal  axis  of  the  aircraft  is  parallel, to  the  airstream  at  all 
times . 

The  root  locations  possible  under  the  adopted  scheme  were  not  exactly 
as  desired  but  represented  the  most  acceptable  in  overall  considerations. 
The  short  period  pair  are  not  well  damped  but  the  damping  ratios  obtainable 
are  fairly  consistent.  For  flight  conditions  1,  2,  5,  and  6  the  frequency 
is  lower  than  that  desired  but  this  was  characteristic  of  these  conditions 
for  all  schemes .  Also  quite  important  is  the  fact  that  no  real  roots  were 
unstable  for  the  range  of  desirable  JC.  with  the  exception  of  condition  3„ 
This  condition  possesses  a  real  root  segment  in  the  right  half  plane ,  but 
short  period  root  location  could  be  effected  with  the  unstable  root  still 
close  to  the  origin.  For  the  range  of  gain  necessary  for  adaptive  control 
the  phugoid  pair  does  not  become  oscillatory  in  any  flight  condition. 
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In  flight  denditions  1  and  5  the  system  exhibits  relative  degrees  of 
insensitivity  to  gain  changes .  This  is  characteristic  of  these  conditions 
mainly  because  of  the  low  q  involved,  while  for  conditions  3  and  U   the 
system  is  very  sensitive  to  gain  changes.  These  are  the  conditions  with 
a  high  q.  However,  N  plus  filter  1,  with  the  mentipned  drawbacks,  still 
represents  the  best  scheme  and  as  such  was  selected  for  the  Adaptive 
Controller  phase  of  the  study. 

The  next  scheme  was  that  of  N  with  filters  2  and  3.  This  scheme  was 
capable  of  providing  a  system  for  self  adaptive  control.  However,  in 
comparison  with  the  N  plus  filter  1  scheme  it  was  less  desirable.  The 
range  of  gain  values  required  was  greater,  and,  even  though  achievable 
damping  ratios  were  better,  the  frequency  of  the  short  period  pair  was 
lower  resulting  in  a  slower  response.  The  phugoid  pair  held  their  positions 
on  the  negative  real  axis  and  then  became  unstable  very  quickly  with  in- 
creasing IC.  .  Thus  the  phugoid  pair  was  oscillatory  for  just  a  very  small 
range  of  gains  and  obviously  at  a  very  low  frequency.  In  fact  this  pair, 
in  general,  were  the  roots  which  provided  the  stability  limit. 

As  in  the  case  of  the  other  investigations,  flight  conditions  3  and 
U  were  the  most  sensitive  to  gain  changes  0  The  gains  could  be  kept  low 
enough  in  each  of  these  conditions  to  prevent  any  spiral  divergence 
condition  from  arising,  and  while  the  damping  on  these  short  period  pair 
was  acceptable  their  companion  frequency  was  too  high. 

The  perturbation  scheme,  using  filter  3»  proved  to  be  even  more 
unsatisfactory.  It  contributed  loci  essentially  the  same  as  those  with 
filter  2  but  they  were  much  more  insensitive  to  gain  changes.  The  response 
was  slower  and  in  general  filter  3  was  no  Improvement  on  filter  2. 

The  next  scheme  used  in  the  analysis  was  N  with  r.  This  scheme 
was  also  acceptable  for  self  adapting  control  and  required  no  compensation 
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filters.  The  damping  ratios  are  good,  the  frequency  is  low,  and  the 
range  of  gains  required  is  large.  However,  the  principal  drawback  might 
be  the  presence  of  aerodynamic  roll  cross-coupling  effects  which  enter 
because  of  the  r  term,  and  make  coordinated  maneuvers  somewhat  difficult. 
Since  only  the  feasibility  of  this  control  was  investigated,  no  attempt 
to  apply  it  to  the  adaptive  controller  was  forthcoming. 

The  investigation  of  this  scheme  does  give  some  information  in  rather 
an  indirect  manner.  While  the  lateral  accelerometer  is  to  be  placed  at 
or  displaced  from  the  aircraft  center  of  gravity,  the  c0g0  may  move.  The 
analysis  of  this  scheme  shows  that  such  movement  affects  principally  the 
frequency  of  the  short  period  pair  and  not  their  damping  ratio.  This  may 
provide  some  limits  on  the  amount  of  e.g.  shift  a  system  using  this  scheme 
can  tolerate . 

The  last  response  analyzed  was  that  of  yaw  rate  feedback.  This  scheme 
gave  good  damping  ratios  and  a  reasonable  range  of  gain  required 0  The 
frequency  of  the  short  period  pair,  while  low,  is  acceptable.  However, 
the  scheme  is  essentially  unacceptable,  with  an  integrated  autopilot, 
from  an  aerodynamic  standpoint .  Tfte  roll  coupling  could  prevent  coordinated 
turns.  As  the  aircraft  banked  into  a  turn  the  yaw  rate  introduced  might 
oppose  the  change  in  heading.  Considering  a  three  degree  of  freedom 
aircraft  this  scheme  is  quite  satisfactory  but  with  a  five  degree  of  freedom 
"  aircraft  it  may  be  unacceptable . 

As  an  integral  portion  of  the  study  the  analog  computer  was  used  to 
verify  the  predicted  performance  obtained  from  the  root  locus  analysis . 
The  three  degree  aircraft  equations  of  motion  as  well  as  all  loop  transfer 
equations  were  mechanized  for  the  analog  computer  as  indicated  in  Appendix 
B.  Transient  characteristics  were  investigated  for  all  lateral  response 
schemes  with  the  exception  of  N  plus  Kir   . 
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The  major  reasons  for  this  portion  of  the  analysis  were  to  determine 
response  characteristics  not  discernible  from  the  root  locus  analysis. 
The  inputs  were  1  deg.  rudder  deflections  (£0  and  1  deg.  p  ,  or  sideslip 
angle,  gusts.  The  loop  gain  was  varied,  as  in  the  root  locus  approach,  to 
establish  the  optimum  gain  for  each  flight  conditions.   In  addition  the 
value  of  gain  at  the  stability  limit  was  obtained.  These  values  agreed 
rather  well  with  those  obtained  in  the  root  locus  study. 

The  transient  characteristics,  as  determined  from  the  analog  study, 
are  tabulated  in  Tables  XIII,  XIV,  and  XV.  Comparing  Table  XI  with  Table 
XIII  for  the  case  of  r  feedback  it  may  be  noted  that  for  flight  condition 
3  the  second  order  approximation  is  quite  good.  However,  for  the  other 
flight  conditions  particularly  damping  ratio  correlation  is  poor.  This  is 
due,  mainly,  to  distortion  of  the  transient  response  due  to  real  root 
influence.  This  fact  can  be  illustrated  by  referring  to  Fig.  4-3  which 
illustrates  the  most  acceptable  r  responses  obtained  at  each  flight  condition. 
Since  the  real  root  distortion  of  the  basic  second  order  response  was  present 
in  every  flight  condition,  the  existence  of  errors  in  determination  of  damping 
ratios  is  probable  0  Therefore ,  the  characteristic  response  is  not  approx- 
imated by  a  second  order  system  for  the  flight  conditions. 

For  the  N  with  filter  2  scheme  the  second  order  approximation  was  not 
totally  free  of  distortion.  The  transient  response  data,  for  the  most 
acceptable  N  response  as  shown  in  Fig.  44,  did  not  conform  to  the  estimated 
values  of  damping  ratio  and  frequency.  Again  this  was  the  result  of  real 
root  influence  on  the  transient  characteristics.  Resultant  distortion  of 
the  oscillatory  response  of  the  short  period  roots  made  data  correlation 
impossible  or  inadequate . 
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TABI£  XIII 

Analog  Computer  Second  Order  Approximations  of  Transient 

Characteristics 

Yaw  Rate  Feedback   (No  Filter) 
Inputs  1  deg  L  or  1  deg  p  Gust 

Flight                 *r                        $                       f(cps)  T(sec) 

Cond.       deg/deg/sec       Est.       Actual       Est.  Actual  Est.       Actual 

1  0.2  0.27  0.244  4.1 
0.4  0.44  0.244  4.1 
0.6  0.17  0.263  3.8 
0.8  0.18  0.25  4.0 
1.0  0.18  0.263  3.8 
1.5              0.62        0.23      0.302  0„23  3.32        4.g 

2.0                               0o25  0o228  4.4 

2  0.5                               0o42  0o278  3.6 

0.75                            0.44  0.286  3.5 

1.0              0.7          0o5        0.286  0.27  3.4B         3.7 

1.25                             0.6  0o27  307 

1.5                                .42  

2,0 0q4  -~- 

3  0.1                              0.47  O0588  1.7 

0.15                            0o31  0o588  1„7 

0.2              0.35         0.42      0.7  0.625  1.43        1.6 

0.25                            0.405  0o588  1.7 

0.3                              0o42           •  0.588  1.7 

0.35                            0.5  0o588  107 
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TABIE  XIII  (Cont.) 

Flight  Kr  J  f(cps)  T(sec) 

Cond.       deg/deg/sec       Est.       Actual       Est.       Actual  Est.        Actual 

U  0.2  0.515       0o6         0.475         0.77  2.1  1.3 

0.5  0.6  

0.8 


0.5 
1.0 
1.5  0.52 


H.F. 
osc<<> 
2  cps 

0.28 

0o27 

3.7 

0o2l 

0.278 

306 

0.28 

0.318 

0^278 

3.K 

3.6 

0-1      Response  too  distorted  to  provide  good  estimate 
on  "$    and  f . 
Damping  comparable  to  1  overshoot  and  low 
undershoot . 
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TABLE  XIV 

Analog  Computer  Second  Order  Approximations  of  Transient 

Characteristics 


N     Plus 
7 

Filter 

2 

Inputs  1  deg  bp 

or  1  deg  p  Gusts 

Flight 
Cond. 

V 

deg/g 

Est. 

Actual 

f 
Est. 

(cps) 
Actual 

T 
Est. 

(sec) 
Actual 

1 

10 

0.32 

| 

"    0.256 

3.9 

20 

0.32 

0.25 

4.0 

30 

0.4-2 

0o25 

4.0 

40 

0.475 

0.244 

4.1 

50 

0.475 

.  0.244 

4.1 

60 

0.54 

0.227 

4.4 

70 

0.54 

0.222 

4.5 

80 

0.7 

0.22 

405 

90 

0.4 

0.7 

0.222 

0.22 

4.5 

405 

2 

1.5 

0.2 

0.286 

3.5 

2.0 

0.2 

0.286 

3.5 

2.5 

0.18 

0.286 

3.5 

3.0 

0.2 

0.294 

1 

5 

3.4 

4.0 

0.195 

0.294 

3.4 

8.0 

0.21 

0.303 

3o3 

12 

0.22 

0.333 

3.0 

15 

0.25 

* 

0.303 

3.3 

20 

0.24 

0.303 

3.3 

23 

0.44 



0.222 



4.5 

— 

103 


TABI£  XIV   (Cont.) 


Flight 
Cond. 

deg/g 

* 

Est. 

Actual 

f 

Est. 

(cps) 
Actual 

T   (sec) 
Est.          Actual 

3 

1.5 

1.75 

2.0 

2.25 

0.23 
0.27 

0.285 
0.285 

0.5 
0.526 
0.541 
0.541 

2.0 
1.9 
1.85 

1.85 

2.5 

0.31 

0.32 

0 

.51 

0.5 

1.96 

2.0 

4 

0.3 
0.5 
0.7 
0.9 

0.185 
0.23 
0.32 
0.34' 

0.345 
0.357 

0.333 
0.313 

2.9 
2.8 
3.0 
3.2 

1.4 

o.u 

0, 

.365 

2.74 

1.5 

Diverges 

5 

15 
20 
25 
30 
35 
40 
45 

0.1 

0.09 

0.14 

0.12 

0.14 

0.12 

0.14 

0.266 

0.25 

0.25 

0.25 

0.25 

0o25 

0.238 

3o9 
4.0 
4„0 
4.0 
4.0 
4.0 
4„2 

50 

0.154 

0.12 

0, 

.27 

0.238 

3.7 

402 

6' 

0.3 
0.5 
0.7 

0.9 
1.1 

0.02 

0.06 

0.02 

0.02 
0.02 

0.27 

0.263 

0.27 

0.27 
0.27 

3.7 

3.8 

3.7 

3.7 
3.7 

1.3 

0.11 

0.02 

0 

.27 

0.278 

3.7 

3o6 
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TABI£  XV 

Analog  Computer  Second  Order  Approximations  of  Transients 

Characteristics 


N  Plus  Filter  1 
Inputs  1  deg  &u  or  1  deg  p  Gust  (High  Freq.  Values  in  Parens) 

Flight     \r  $  f  (cps)        T  (sec) 

Condi    deg/g    Est.   \ctual    Est.   \ctual   Fst.   Actual 


1 

10 

0.23 

0.25 

4.0 

15 

Q.27 

0.278 

3.6 

20 

0.312       0.37        0.288         0.263         3.48         3.8 

25 

0.5 

0.286 
(1.67) 

3.5 
(0.6) 

30 

0.5 

0.2 
(1.43) 

5.0 

(0.7) 

37 

0.6 
(0.3) 

0.278 
(1.43) 

3.6 

(0.7) 

40 

0.6 
(0.3) 

0.333 
(1.43) 

3.0 
(0.7) 

45 

— — 

— .— 

- 

(0.2) 

(1.33) 

(0.75) 

50 

(O.H) 

(1.43) 

(0/7) 

55 

___ 

— — 

■■i  .i 

(0.1) 

(1.33) 

(0.75) 

60 

(0.085) 

(lo33) 

(0.75) 

65 

-—. 

... 

■■  - 

(0o02) 

(1025) 

(0.8) 
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TABIE  XV   (Cont.) 


Flight 
Cond  . 

deg/g 

Est. 

Actual 

f 

Est. 

(cps) 

Actual 

T   ( 
Est. 

sec) 
Actual 

2 

1.6 

0.17 

0.286 

3o5 

2.0 

0.185 

0.286 

3.5 

205 

0.18 

0.286 

3.5 

3.0 

0.21 

-    0.29 

3045 

3.5 

0.215 

0.294 

. 

3.4 

4.0 

0.21 

0.303 

3o3 

4.5 

0.22 

0o294 

3.4 

540 

0.45 

Q.24 

0o35 

0.294 

2085 

3.4 

5.5 

0.24 

0.294 

3.4 

3 

0.5 

0.22 

0.265 

0.73 

0.625 

1.37 

1.6 

0.75 

0.33 

0.77 

lo3 

1.0 

0.35 

0.91 

1.1 

1.25 

0.225 

1.11 

0.9 

1.5 

0.0 

1.25 

0.8 

4 

0.3 

0.275 

0.385 

2.6 

0.5 

0.295 

0.435 

2.3 

0.6 

0.4 

0.455 

2.25 

0.7 

0.29 

0.4 

0.605 

0o476 

1.65 

2.1 

o.e 

0.4 

0.5 

2.0 

0.9 

0.47 

TJ.526 

1.9 

1.0 

0.62 

0.^66 

2.15 

1.72 

0.0 

1.22 

0.82 

106 


TABI£  XV   (Cont.) 


Flight 
Cond. 

V 

deg/g 
10 

Est. 

Actual 

f 
Est. 

(cps) 
Actual 

0.278 

T   ( 

Est. 

sec) 
Actual 

5 

0.13 

3.6 

15 

0.14 

0.278 

3.6 

20 

0.25 

0.14 

0.32 

0.278 

3.12 

3.6 

25 

0.17 

0.286 

3o* 

,   30 

0.21 

0.286 

3.5 

35 

0.3 

0.278 
(1.33) 

306 
(0.75) 

40 

0*3 

(1.33) 

(0o75) 

45 

(0.05) 

(1.33) 

(0.75) 

6 

2.0 

0.17 

0.308 

3.25 

3.0 

0„185 

0.333 

3.0 

4.0 

0.24 

0o333 

3.0 

5.0 

0.42 

0.29 

0.382 

0.345 

2.62 

209 

6.0 

0.29 

0.364 

2.75 

7.0 

0.4 

0.4 

2.5 

8.0 

0.5 
(0.15) 

(1.33) 

(0.75) 

9.0 

(0.05) 

(1.25) 

(0.8) 

10  7 


r  (deg/sec) 

-3.0 

0  

3.0 

-20.0 


20.0 

•180.0 
0  


180.0 


-20.0 


0 


20.0 


-15.0 


15.0 
■25.0 
0  ~ 


25.0 


Lime  Scale:  1  div-1  sec 
Fig.  43 


Cond.  1 


Cond.  2 


Cond.  3 


Cond.  4 


Cond.  5 


Cond.  6 


£r=1  deg(+) 


.•dr  frame  Response,  r,  to  Command  Judder  Input  at  Optimum  Loop 
Gain  for  Each  flight  Condition. 
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Ny  (g) 

0.004 
0  


-  0.004 


0.02 


0 


-0.02 


0.1 


0 


-  0.1 


0.2 


0.2 
0.01 


0 


-0.01 


0.1 


0.1 


/  ,u  I  I  I  I  I 


\  \  \  M  \ 


Time  Scale:  1  div  =1  sec 


Fig.  hh 


Jondo  1 


Cond.  2 


^ond.  3 


Cond.  4 


Cond .  5 


Cond.  6 


£r-1  deg(V) 


Airframe  itesponse,  N«,  to  Command  Judder  Input  with  Filter  2 
Compensation  and  at  Optimum  Loop  Gain  for  .^ach  Flight  Condition. 
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In  the  final  scheme,  N  with  filter  1,  studied  on  the  analog  computer 
the  N  response  obtained  was  closer  to  that  predicted  from  the  root  locus. 
While  the  second  order  approximation  was  not  exact  it  was  quite  close  for 
all  six  flight  conditions.  This  may  be  noted  from  a  comparison  of  the  data 
in  Tables  VI  and  XV0  Fig,  45  contains  the  characteristic  response  of  N  , 
at  the  most  acceptable  values  of  K„  ,  for  each  flight  condition.   In 
general,  the  damping  was  less  and  the  frequency  lower  than  the  second  order 
^root  locus  predictions.  However,  the  higher  frequency  oscillation  discussed 
in  Chapter  VII  does  make  itself  evident 0  This  data  is  included  in  Table 
VII 0  The  lower  damping  ratio  and  frequency  represent  the  contributions 
from  roots  in  regions  adjacent *to  the  origin. 

It  must  be  concluded  that  the  analog  study  definitely  verified  the 
choice  of  the  N  with  filter  1  scheme  for  adaptive  control.  The  oscil- 
latory nature  of  the  response  essentially  permitted  a  more  rapid  control 
of  damping  ratio  as  a  function  of  variable  loop  gain0  The  real  root 
distortion  of  the  dominant  complex  roots  in  the  case  of  filter  2  and  in  the 
use  of  r  as  a  lateral  response  were  detrimental  in  provision  of  a  rapid 
indication  of  varying  ^P  0 

Furthermore,  this  analysis  demonstrated  the  feasibility  of  Adaptive 
Control  of  this  system  over  its  flight  regime,  using  gain  changes  only, 
and  indicated  that  this  was  possible  with  several  schemes  of  feedback. 
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Ny(g) 

0.004 


-0.004 
0.02 
0  

-  0.02 

0.1 

0  

-0.1 

0.2 


0 


0.2 

0.01 
0  — 


0.01 

0.04 
0  — 


-0.04 


ConcL  1 


Condo  2 


Cond.  3 


Cond.  4 


Cond.  5 


Cond.  6 


Time  Scale:  1  div.=  1  sec 


Fig.  45 


Sr= 


R=l  deg(V) 


Airframe  Hesponse,  Ny,  to  Command  Judder  Input  with  Filter  1 
Compensation  and  at  Optimum  Loop  Gain  for  x-ach  I light  Condition. 
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Chapter  XII 
Analysis  of  the  Self  Adaptive  Schema 

On  the  basis  of  the  stability  analysis  conducted  on  the  aircraft, 
the  lateral  airframe  response  corresponding  to  N  with  lead  filter 
compensation,  filter  1,  was  selected  as  the  response  to  be  controlled. 
Values  of  IL,  were  chosen  as  those  variable  gain  values  which  permitted 
the  criterion  value  to  minimize .  In  terms  of  damping  ratios ,  this 
corresponds  to  a  f  equal  to  0o7  or  if  not  obtainable,^  equal  to  a  max- 
imum value  0  Table  XVI  contains  the  values  of  K~  optimum  and  J*  corres- 
ponding to  IL.  optimum  determined  by  results  of  the  analog  computer  stability 
study 0 

Input  disturbance  to  the  airframe  was  supplied  by  a  short  duration 
pulse  input  inserted  by  signal  generator  into  the  loop  at  the  input  side 
of  the  variable  gain  pot0  The  pulse  was  below  the  pilot's  threshold  value 
of  0o05  g  lateral  acceleration. 

Results  of  the  adaptive  scheme  as  applied  to  the  selected  lateral 
response  are  contained  in  Table  XVII.  These  data  are  basically  arranged 
with  regard  to  flight  condition  and  may  be  discussed  in  this  manner. 

The  landing  condition,  condition  1,  was  tested  and  found  to  be  extremely 
sluggish  in  adaptive  corrective  response „     While  the  optimum  gain  selected 
from  the  analysis  should  have  been  JL.    equal  to  25,0  deg/g,  the  system 
barely  corrected  from  the  initial  1L.    of  0.5  in  a  period  of  two  minutes. 
More  success  was  obtained  by  letting  IL-  equal  10  initially.  In  this  case, 
the  adaptation  was  more  satisfactory  but  a  1C.  equal  to  30.0  deg/g  was 
selected  by  the  scheme.  Finally,  an  initial  JL.  of  4-0.0  deg/g  was  incor- 
porated and  adaptation  served  to  produce  no  change  in  the  variable  loop 
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TABLE  XVI 
Values  of  1L.    for  Desired  Optimum  and  Corresponding  j 


\ 


Flight 
Cond. 

V 

optimum)  deg/g 

f 

1 

25.0 

0.5 

2 

5.5 

0.24 

3  . 

1.0 

0.35 

4 

* 

0.9 

0.47 

5 

30.0 

0.21 

6 

7.0 

0.4 
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gain.  This  testing  was  carried  out  with  zero  reference  signal.  Testing 
at  ^  0.4-  volt  reference  signal  indicated  that  the  flight  condition  was 
essentially  insensitive  to  the  polerity  of  the  reference  but  did  tend  to 
adapt  very  slowly  for  the  ♦  0oi  volt  reference.   In  all  cases  tested,  the 
maximum  lateral  acceleration  was  well  within  the  threshold  of  the  pilot. 
The  system  recovered  when  purposely  given  an  unstable  loop  gain  as  an 
initial  value .  The  final  K~,  sought  by  the  mechanism  was  varied  about 
the  optimum  value  selected  from  the  root  locus  plots . 

Plight  condition  2  exhibited  satisfactory  tendencies  while  using  a 
♦  0.4.  volt  reference  signal.  Adaptation  for  zero  or  negative  referencing 
was  either  prohibitively  slow*or  unsatisfactory.  Under  the  positive  ref- 
erence signal,  the  IC.  sought  by  the  system  was  at  all  times  between  the 
values  of  3o0  to  6o0  deg/g.  This  compares  with  an  optimum  JL.    of  5.5 
deg/go  At  K~  equal  to  0o5  initially,  adaptation  was  slow.  The  system  was 
able  to  seek  the  optimum  gain  region  with  conditions  of  stable  and  unstable 
loop  gains  initially  imposed.  An  investigation  of  the  behavior  under 
increased  pulse  disturbance  resulted  in  the  system  seeking  the  3.0  value. 
The  maximum  lateral  acceleration  experienced  in  the  stable  testing  above 
was  0.02  g. 

Flight  condition  *3  as  with  the  previous  flight  conditions  performed 
most  satisfactorily  with  a  positive  reference  value.  With  negative  reference 

applied,  the  loop  gain  correction  tended  to  be  applied  in  the  wrong  direction 
while ,  with  a  zero  reference ,  corrective  action  approached  that  of  the 
positive  reference.  The  maximum  stable  N  experienced  was  0.02  g  with  an 
increased  disturbance  amplitude.  The  average  N  was  0.01  g.  Since  this 
flight  condition  was  the  only  one  which  was  divergently  unstable e  it 
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furnishes  the  only  information  on  stable  recovery  from  this  condition. 
The  system  was  capable  of  adjusting  from  divergent  instability  to  a 
IL.  equal  to  lo0  deg/g  in  5  sec.  With  the  positive  reference  signal, 
the  system  sought  to  adapt  the  variable  loop  gain  to  a  value  of  IL,  from 
0.5  to  0o6  deg/g. 

The  most  sensitive  flight  condition,  condition  4.,  adapted  quickly, 
usually  after  a  maximum  of  two  optimizing  periods.  In  this  case  the 
maximum  N  observed  during  adaption  was  the  pilot  threshold,  0o05  g. 
Both  positive  and  zero  reference  signals  provided  the  most  satisfactory 
response.  The  system,  with  K~  initially  0o5  deg/g,  tended  to  remain  at 
this  value.  After  an  excessive  amount  of  time  the  selected  gain  sought 
to  increase  to  a  greater  value .  This  caused  a  decrease  in  f   and  the 
resultant  immediate  adaptation  to  a  Kj.  of  0.5  to  0o9  deg/g  as  compared 
to  the  optimum  of  0.9  deg/g 0  Since  this  lower  end  value  also  constitutes 
the  residual  loop  variable  gain  there  is  no  data  to  support  the  tendency 
of  the  adaptive  scheme  to  remain  at  this  value  rather  than  decrease 
further.  As  a  result  of  the  airframe  sensitivity  at  this  flight  condition, 
there  was  a  noticeable  tendency  for  the  system  to  attempt  to  hunt  or  drift 
about  values  of  IL.  to  which  it  adapted  0  This  tendency  was  actually  noticed 
to  produce  a  more  positive  adapting  action  since  the  system  tried  to 
continuously  correct  rather  than  sluggishly  allow  considerable  error  to 
accumulate  0 

Flight  condition  5»  a  sluggish  flight  condition,  exhibited  satisfactory 
tendencies  with  positive  or  zero  reference  signals  applied.  Divergent  or 
unstable  tendencies  were  not  present.  The  system  was  capable  of  recovering 
from  near  unstable  or  very  low  $   and  high  frequency  conditions .  Since  the 
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maximum  gain  possible  with  the  variable  pot  was  4.0.5  deg/g  it  was  not 
possible  to  drive  the  aircraft  unstable  and  inspect  recovery  tendencies . 
For  IL.  initially  equal  to  0.5  deg/g  the  system  should  have  adapted  to  the 
optimum  value  of  30.0  deg/g.  Due  to  the  relative  insensitive  aircraft 
response  to  the  input  disturbance  however,  very  little  corrective  tendency 
was  noted.  After  a  period  of  7  minutes  the  variable  gain  had  adjusted 
from  0.5  to  only  0.7  deg/g0  An  initial  fL.  of  20  o0,  which  coincided 
with  the  optimum  gain  value,  was  imposed  upon  the  system.  The  controller 
did  not  attempt  to  adapt  from  this  position.  Data  is  not  conclusive 
enough  to  conclude  that  the  scheme  sensed  the  optimum  condition  since  in 
the  sluggish  condition  of  flight  the  value  of  J3  varies  imperceptibly  with 
a  small  variance  in  JL.    .  Under  the  stable  conditions  of  test  and  maximum 
N  achieved  was  0.02g. 

Finally,  condition  6  with  an  optimum  K~  of  5.0  deg/g  exhibited 
fairly  satisfactory  tendencies  at  reference  signals  equal  to  zero  and 
positive  0o4.  volts.  However,  the  K-.  sought  was  in  the  region  1.0  to  2o0 
deg/g.  With  an  initial  1L,    of  3o0  the  tendency  of  the  system  was  to  select 
a  variable  gain  value  of  6.0  deg/g  0  Increasing  the  disturbance  magnitude 
with  the  same  initial  conditions  tended  to  cause  a  selection  of  variable 
gain  in  the  region  of  3.0  deg/g  with  hunting  about  this  value.  Imposing 
an  initial  condition  of  5.0  deg/g  with  the  increased  disturbance  r«iulted 
in  hunting  about  the  value  K-  equal  to  6.0  deg/g 0  The  maximum  fty  achieved 
was  no  greater  than  0.02g  during  this  testing.  As  In  condition  U   tests,  it 
was  noticed  that  adaption  was  relatively  faster  since  airframe  response  was 
generally  more  sensitive  to  a  disturbance . 
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J  2  sec. 
error   dt  for 
o 
1  deg.  bu  versus  variable  gain  JL.     ,  for  each  flight  condition  tested 

where  | error  |  equals  |N  plus  filter  1  j  and  is  measured  at  the  output 
of  filter  1  to  take  advantage  of  a  larger  signal 0  Comparison  of  these 
figures  to  Fig.  38  indicate  that  the  second  order  variance  was  not 
reasonably  attained.  The  most  resemblance  was  attained  in  those  cases 
where  the  real  root  contributions  were  less  distorting  to  the  short 
period  characteristics .  It  may  be  noted  that  miniraums  occur  in  each  error 
function  in  the  general  vicinity  of  the  predicted  optimum  IL,  .  Further- 
more ,  several  of  the  error  integrals  display  minimums  at  values  other 
than  at  the  predicted  optimum  JL.     indicated  in  Table  XIII0  Fig.  46 
indicates  a  minimum  in  the  vicinity  of  IC.  equal  to  22.5  deg/g  while  there 
are  also  exhibited  M flats1'  or  negligible  change  in  the  error  integral  on 
either  side  of  the  minimum.  Fig.  47  for  flight  condition  2  indicates  a 
minimum  at  about  K_  equal  to  5o0  deg/g.  This  coincides  with  the  predicted 
minimum.  Condition  3  shown  in  Fig.  48  indicates  a  definite  minimum  at 
JL.    equal  to  1025  deg/g  but  a  region  of  maximum  error  at  K~  equal  to 
0.7  deg/g  the  predicted  optimum  value.  For  all  intents  and  purposes  the 
function  exhibits  a  minimum  at  IL.  equal  to  0o5  deg/g,  the  residual  pot 
gain,  should  the  system  attempt  to  initially  adapt  from  this  value.  Fig. 
49,  the  most  critical  flight  condition,  exhibited  a  satisfactory  variance 
of  the  error  function  with  varying  loop  gain0  The  curve,  at  all  gain 
values ,  slopes  so  that  increases  or  decreases  of  gain  from  the  minimum 
serve  to  always  increase  the  criterion  value.  Fig.  50,  condition  5. 
indicates  that  the  error  exhibits  flat  regions  from  0o5  to  35.0  deg/g. 
This  flat  region  coincides  with  the  minimum  values  of  the  integral.  The 
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final  flight  condition  indicated  in  Fig.  51  exhibited  an  essentially  flat 
region  from  the  residual  gain  K-.  equal  to  0.5  deg/g  to  approximately 
Kj,  equal  to  2.0  deg/g.  The  error  integral  then  achieved  a  maximum  and 
thereafter  a  minimum  value  of  K„  equal  to  6,0  deg/g.  Again,  in  this 
flight  condition,  there  were,  as  far  as  the  corrective  logic  could 
discern,  two  seekable  minimums  which  depended  upon  the  initial  value  of 
K«  imposed,  or,  the  value  of  variable  gain  with  relation  to  the  criterion 
maximum  mentioned  above , 

Therefore,  in  the  case  of  flight  condition  1,  the  landing  condition 
and  one  in  which  maximum  stability  is  desired  from  a  control  standpoint, 
it  was  ascertained  that  adaptation  was  sluggish  and  slow  due  to  relatively 
small  values  of  the  error  criterion  as  a  result  of  the  input  disturbance . 
Furthermore ,  the  flat  regions  of  criterion  value  caused  little  or  no 
variation  with  gain  and  therefore  the  criterion  was  actually  insensitive 
to  gain  in  these  regions .  These  tendencies  will  result  in  the  adaptive 
controller  remaining  in  the  region  or  perhaps  slightly  drifting  in  variable 
gain  value  until  either  an  increase  or  decrease  in  criterion  value  causes 
further  adjustment 0  Admittedly,  the  change  is  oscillating  damping  ratio 
between  these  regions  is  small  and  this  may  be  noted  from  the  root  locus 
plot  in  Figc  6.  Since,  as  seen  from  the  figure,  the  movement  of  the  short 
period  root  from  the  complex  pole  is  slight  over  the  gain  region  0,5  deg/g 
to  4-5  deg/g. 

Therefore,  it  may  be  concluded  that  the  ability  of  the  system  to  seek 
the  flat  or  minimum  regions  of  the  criterion  was  considered  adequate  while 
its  ability  to  seek  and  remain  at  the  minimum  value  of  the  criterion  was 
inadequate  0  Furthermore ,  it  may  be  concluded  that  relatively  small  values 
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of  criterion  signal  resulted  in  a  small  adaptive  error  or  drive  signal 
to  the  servo  gain  changes  which  made  corrective  action  sluggish  even  at 
high  servo  gains .  The  positive  reference  signal  was  adequate  for 
adaptation  and  maximum  Ny  was  acceptable.  Finally,  the  system  exhibited 
the  ability  to  recover  from  oscillatory  instability  and  adapt  to  an  accept- 
able response . 

Flight  condition  2  effectively  adapted  in  the  regions  of  K„  equal 
to  3.0  to  6.0  deg/g.  This  agrees  with  the  region  of  minimum  values  of 
the  error  criterion.  Due  to  the  flat  region  exhibited  by  the  low  gain 
region  of  the  error  criterion,  adaptation  from  the  residual  1L.     of  0o5 
deg/g  was  unacceptably  slow.  "This  can  be  attributed  to  the  fact  that  the 
adaptive  error  or  criterion  value  remained  essentially  constant.  The  slope  of 
the  criterion  curve  for  values  of  IL.  greater  than  the  optimum  value  of  5.0 
deg/g  was  adequate  to  determine  quickly  an  increase  in  the  criterion  and 
apply  the  necessary  correction.  The  performance  under  increased  pulse 
disturbance  was  satisfactory  and  desirable  under  the  existing  conditions 
since  the  increased  pulse  would  supply  more  positive  directional  information 
at  the  gain  region,  IL.  equal  to  0.5  to  2Q5  deg/g.  The  system  was  also 
capable  of  recovering  from  unstable  initial  conditions  and  did  so  adequately. 
The  maximum  Ny  achieved  during  stable  adaptation  was  acceptable  and  well 
below  the  pilot's  sensing  threshold. 

Therefore,  it  may  be  concluded,  that  for  flight  condition  2  the 
positive  reference  signal  was  satisfactory;  adaptation  at  all  initial  values 
of  Kj.  greater  than  or  equal  to  3.0  deg/g  was  satisfactory;  and,  maximum 
airframe  stresses  were  tolerable.  Furthermore,  it  may  be  stated  that 
performance  of  the  system  at  initial  values  of  K„  less  than  3o0  deg/g  was 
slow  and  essentially  ineffective . 
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In  the  case  of  flight  condition  3  the  positive  reference  signal  of 
0«X  was  again  found  to  be  the  most  effective  and  acceptable.  The  system 
showed  the  ability  to  recover  adequately  from  divergent  instability  in 
approximately  one  optimizing  interval.  This  particular  flight  condition 
indicated  compatibility  with  the  maximum  and  minimum  values  of  the  error 
criterion.  When  forced  to  adapt  from  an  initial  JL.     of  0.5  deg/g  the 
system  sensed  a  minimum  and  was  unable  to  attain  the  desirable  minimum 
at  approximately  a  IL.  of  1.2  deg/g.  On  the  other  hand,  when  recovering 
from  divergent  instability,  the  system  sought  to  minimize  in  the  preferred 
location.  In  the  other  requirements  the  flight  condition  appeared  to  be 
adequately  handled  by  the  adaptive  logic. 

It  may  be  concluded  then  that  the  system  was  satisfactory  in  adapting 
speed,  recovery  from  divergent  instability  and  operation  with  positive 
reference  signal.  The  data  also  conclusively  shows  that  when  approaching 
the  flight  condition  from  higher  gains  than  IL.  equal  to  1.2  deg/g  the 
system  will  adequately  adapt  to  this  minimum. 

Flight  condition  h.   contains  the  most  satisfactory  performance  data 
with  regard  to  operation  of  the  adaptive  controller.  The  system  sought 
a  minimum  between  the  values  of  K~  equal  to  0.5  deg/g  and  0.9  deg/g.  The 
error  criterion  variance  with  gain  indicates  that  the  minimum  Ny  shoald 
occur  at  a  IL.  of  0.9  deg/g.  The  minimum,  however,  is  critical  since  the 
value  of  the  error  rapidly  increases  on  either  side.  The  hunting  between 
the  values  stated  above  can  be  most  intelligently  reasoned  by  noting  that 
the  adaptive  error  would  be  large  and  the  necessary  gain  variance  is  small e 
Slight  overdriving  of  the  adaptive  positioner  could  result  in  the  system 
sensing  an  error  increase  and  overdriving  the  minimum  which  again  results 
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in  an  error  increase  which  results  in  overdriving  the  minimum  in  the 
opposite  direction.  Admittedly  non-linear  variance  of  the  potentiometer 
at  these  lower  gain  values  increases  the  system  sensitivity  but  apparently 
not  enough  to  allow  the  system  to  remain  at  the  minimum  error  position. 
As  previously  stated,  this  action  was  not  detrimental  to  the  stability 
of  the  aircraft  but  served  to  eliminate  the  sluggish  nature  of  response 
encountered  in  conditions  1  and  2. 

It  may  be  concluded  that  with  regard  to  condition  A   the  adaptive 
system  operated  satisfactorily  with  regard  to  corrective  tendencies .  The 
system  never  was  capable  of  settling  to  the  minimum  value  of  the  error 
criterion  due  to  large  drive  voltages  to  the  positioning  servo  causing 
overshoot  not  completely  compensated  for.  The  system  was  obviously  capable 
*  of  recovering  from  situations  approaching  instability  as  evidenced  by  its 
tendency  to  never  exceed  a  K„  of  approximately  0.9  deg/g.  The  value  of 
the  criterion  increased  so  rapidly  beyond  this  point  that  corrective  action 
was  always  in  the  stabilizing  direction.  Conclusions  about  recovery  from 
unstable  conditions  at  this  flight  condition  are  meaningless  since  the 
aircraft  could  not  possibly  withstand  the  structural  stresses. 

Condition  5  performed  adequately  with  the  positive  reference  signal 
but  conformed  to  the  relatively  insensitive  variance  of  the  error  criterion 
with  variable  gain  by  adapting  in  a  sluggish  manner  at  all  imposed  initial 
conditions.  Since  the  variable  gain  has  very  little  effect  upon  the  detect- 
able response  in  the  region  of  t-   from  0.5  deg/g  to  4.0.0  deg/g,  the 
location  of  the  gain  within  this  region  is  significantly  unimportant. 
The  major  requirement  in  this  instance  is  that  the  departure  of  the  aircraft 
from  this  flight  condition  create  suitable  changes  in  the  error  criterion 
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variance  with  gain  to  permit  selection  of  corrective  information  for  the 
flight  condition  being  approached  which  may,  and  probably  does,  require 
more  exacting  placement  or  positioning  of  the  gain  pot  controlling  JL,  . 

Therefore,  it  may  be  concluded  that  the  tests  at  flight  condition  5 
are  inconclusive  since  practically  any  gain  within  the  region  of  the 
adjustable  gain  pot  results  in  an  undetectable  variation  in  the  error 
criteridn. 

Finally,  in  flight  condition  6,  the  positive  reference  signal  proved 
acceptable  with  regard  to  performance  tendencies .  The  variance  of  the 
error  criterion  with  variable  IL.  exhibited  two  minimum  values  separated 
by  a  maximum.  The  effect  of  this  maximum  was  undesirable  since  it 
conceivably  could  result  in  performance  at  the  flight  condition  at  some 
non-optimum  point  determined  by  the  minimum  sought  at  values  of  low  IL,  . 
This  is  borne  out  by  the  actual  results  of  testing  which  indicated  the 
tendency  of  the  controller  to  minimize  at  two  separate  values  of  gain. 
These  IL.  values  coincide  approximately  with  the  observed  minimums  in 
Fig.  51. 

Therefore,  it  may  be  concluded  that  adaptation  of  flight  condition  6, 
while  being  satisfactory  with  respect  to  adapting  speed  and  positive  ref- 
erence signal,  was  unsatisfactory  in  being  able  to  seek  the  minimum  value 
of  the  error  criterion  under  all  imposed  initial  conditions .  The  controller 
did  conclusively  prove,  however,  that  it  was  capable  of  seeking  a  minimum 
which,  though  not  the  minimum  value  of  the  error  criterion,  was  the  only 
minimum  the  system  could  discern. 

In  summary,  the  analysis  of  the  adaptive  controller  unit  must  be 
considered  to  establish  the  fact  that  the  basic  logic  and  principles  of 
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controller  design  are  satisfactory.  It  is  concluded  that  the  designed 
sequence  and  logic  of  the  controller  is  acceptable  for  application  to  any 
control  process  whose  parameter  or  parameters  may  be  altered  as  the  direct 
result  of  the  minimization  of  some  measurable  variable.  Since  the  adapt- 
ive controller  is  definitely  limited  by  the  necessary  time  sequencing,  the 
period  of  control  that  must  be  exhibited  by  the  adaptive  scheme  is  of 
considerable  importance  and  must  be  studied  and  proved  acceptable  before 
application.  There  was  insufficient  mechanization  available  on  the 
computer  to  enable  test  of  the  adaptive  system  over  some  characteristic 
flight  profile  to  establish  that  the  controller  was  capable  of  continuously 
modifying  the  loop  gain  to  insure  the  most  acceptable,  minimum  error  res- 
ponse. Study  of  the  variation  in  the  value  of  the  absolute  integral  of 
error  indicated  that  the  system  would  most  probably  exhibit  no  difficulty 
In  adapting  from  large  gain  values  to  the  smaller  values  of  K~  required 
for  the  more  sensitive  flight  conditions.  However,  these  same  curves 
indicated  that  in  some  instances  of  low  gain,  false  minima  occurred  because 
of  real  root  incluence  and  these  may  result  in  the  adaptive  system  being 
subjected  to  minimum  values  which  were  stable  but  certainly  did  not  produce 
the  most  acceptable  response. 

Therefore ,  it  is  considered  that  there  is  inadequate  data  to  enable 
prediction  of  whether  or  not  the  adaptive  system  will  successfully  handle^ 
a  particular  flight  profile. 

The  information  obtained  on  the  behavior  of  the  error  criterion  during 
the  variation  of  gain  for  each  particular  flight  condition  indicates  that 
for  the  aircraft  the  selection  of  the  integral  of  the  absolute  error  was 
unsatisfactory.  The  principal  reasons  for  its  inadequacy  were  the  fact  that 
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in  several  Instances  gain  variations  produced  little  or  no  effect  upon 
the  value  of  the  criterion  and  the  fact  that  false  minima  occurred  in 
several  flight  cases .  The  former  was  attributed  to  the  slow  movement 
of  the  roots  from  the  short  period  poles  causing  very  little  change  in 
the  value  of  the  resulting  error.  The  latter  inadequacy  was  directly 
attributable  to  the  influence  of  the  real  roots  over  the  complex  dominant 
roots  at  low  values  of  IL.  . 

Therefore ,  it  is  concluded  that  the  error  criterion  was  inadequate 
in  the  majority  of  flight  conditions  for  the  reasons  given  above .  Since 
the  adaptive  controller  depends  entirely  upon  the  establishment  of  adequate, 
constantly  minimizing  error  criteria,  it  must  be  furthermore  concluded  that 
establishment  of  the  system  for  the  particular  aircraft  tested  was  not 
achieved. 

Finally,  it  must  be  concluded  that  the  investigation  of  all  possible 
methods  of  compensation  and  establishment  of  appropriate  error  criteria 
which  do  properly  behave  has  by  no  means  been  exhausted .  Therefore ,  it 
is  recommended  that  additional  testing  and  analysis  in  these  particular 
areas  be  conducted.  Detailed  recommendations  with  these  thoughts  in  mind 
are  included  in  the  following  chapter. 
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CHAPTER  XIII 
Conclusions  and  Recommendations 
The  stability  analysis  conducted  using  the  various  suggested  lateral 
responses  as  methods  of  providing  a  means  for  stable  adaptive  control 
resulted  in  the  following  conclusions : 

1.  Yaw  rate  and  N  with  filters  2  and  3  were  undesirable  as  control 

y 

variables  since  the  second  order  dominance  of  the  short  period 

pair  of  roots  were  greatly  distorted  by  the  real  roots .  Due  to 

this  fact,  the  system  tended  to  be  less  oscillatory  and  there- 
i 

fore  less  sensitive  to  a  variation  in  &   as  a  result  of  varying 
loop  gain. 

2.  The  aerodynamic  inadequacies  of  yaw  rate  as  a  response  precluded 
its  use  in  an  integrated  control  loop  of  a  flight  control  system. 

3.  The  feasibility  of  the  combined  response  N  with  r  was  established 

Jr 

with  regard  to  provision  of  a  system  whose  characteristics  more 
closely  resembled  those  desired  for  adaptive  control  by  the 
methods  contained  in  the  report. 
'  4-.  N  with  filter  1  compensation  was  the  most  desirous  of  the 

responses  analyzed,  however,  in  several  insensitive  flight 
conditions  variation  of  $   with  varying  JL.    was  inadequate. 
With  regard  to  the  analysis,  design,  and  performance  of  the  adaptive 
controller  the  following  conclusions  are  advanced: 

'10  Operation  and  performance  of  the  logia  and  sequencing  of  the 

adaptive  controller  was  satisfactory  and  design  was  adequate. 
2,  The  error  criterion,   1  |  E  |  dt  ,  chosen  as  the  vehicle  by  which 
the  controlled  response  was  to  be  minimized  was  inadequate  in  many 
flight  conditions  since  it  failed  to  produce  directional,  corrective 
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information  to  the  controller.  In  the  case  of  flight 
conditions  3  and  6  it  exhibited  minimum  values  which  were 
in  addition  to  the  desired  minimum  occurring  at  the  optimum 
response . 
30  Performance  indicated  that  provision  of  the  proper  minimizing 
criterion,  as  might  be  indicated  by  application  of  the  tested 
error  criterion  to  a  second  order  servomechanism,  would  result 
in  adaptation  to  the  specified  minimum  criterion  value . 
The  recommendations  shall  be  primarily  concerned  with  the  further 
analysis  of  the  lateral  loop  response  and  refinements  or  variations  of  the 
Self  Adaptive  Controller.  The  ramifications  of  some  of  the  problems  con- 
cerned with  a  five  degree  of  airframe  will  also  be  pointed  out  as  paths 
for  further  study  to  proceed  upon. 

Certainly  the  study  of  additional  compensation  devices  should  be 
pursued.  The  effect  on  the  root  loci  family  for  N  with  perturbations 
of  the  lead  filter  pole  and  zero  should  be  investigated.  The  use  of  more 
complex  compensation  to  place  complex  conjugate  zeros  at  a  location  to  aid 
in  forcing  the  short  period  pair  of  roots  into  a  desirable  domain  by  high 
gains  may  be  possible  and  worthy  of  study.  Study  should  be  conducted  upon 
the  effect  of  lead  compensation  on  the  N  plus  f  scheme . 

Concerning  the  adaptive  controller,  the  greatest  area  for  study  is  in 
the  selection  of  a  suitable  error  criterion.  A  criterion  must  have  a 
distinguishable  minimum  at  the  optimum  value .  Pref errably  it  should  be 
concave  in  shape  with  a  single,  well  defined  minimum.  To  be  able  to  establish 
this  characteristic  for  an  error  criterion  for  any  complex  servo  would  make 
this  Adaptive  Controller  very  effective*.  The  time  sequence  of  the  controller 
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should  be  studied  further.  The  relative  length  of  the  operating  portions 
of  the  optimizing  period  as  well  as  the  overall  length  of  the  period  need 
investigation  to  choose  their  optimum.  Studies  to  develop  systems  to 
establish  procedures  for  choosing  of  optimum  times  would  be  of  value. 

Further  study  should  be  made  of  the  methods  available  for  introducing 
disturbance  signals  to  the  airframe.  Studies  should  provide  for  invest- 
igation of  the  most  suitable  type  of  disturbance 0  Certainly  the  use  of 
random  noise  to  actuate  the  controller  instead  of  pulsing  the  rudder 
would  be  of  an  advantage.  However,  the  effect  of  this  would  have  to  be 
determined  with  regard  to  adequacy,  magnitude  and  sensitivity. 

Another  area  for  study  would  be  the  gain  changing  servo  loop  itself. 
The  use  of  rate  limiting  on  this  loop  was  proposed  to  enable  it  to  adapt 
faster  yet  be  less  sensitive  to  a  gusting  environment.  Further  study  is 
necessary  to  pick  the  best  value  for  the  threshold  which  changes  the  gain 
of  this  loop. 

Considering  a  five  degree  of  freedom  airframe  many  problems  arise. 
The  solution  of  these  problems,  or  at  least  a  feel  for  their  effect,  would 
be  necessary  before  design  is  considered  final.  The  effect  of  b.  inputs 
must  be  further  investigated.  The  analysis  and  testing  of  the  effect  of 
roll  cross  coupling  would  be  of  great  interest.  The  effect  of  this  coupling 
would  be  of  great  interest.  The  effect  of  this  coupling  when  using  schemes 
like  N  plus  r  should  be  analyzed.  The  effect  of  residual  vibrations  or 
environmental  vibrations  on  instrumentation  should  be  investigated 0  The 
analysis  of  the  airframe  considered  as  flexible,  with  the  consequent 
involvement  of  bending  modes,  would  complete  any  study  of  this  control  loop. 
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APPENDIX  A 
Computation  ffethods  for  Determination  of 
Root  Locus  Plots 

In  order  to  conveniently  compute  the  numerous  root  locus  plots 
required  during  the  analysis,  the  stability  matrix  associated  with 
the  control  loop,  or  its  characteristic  determinant,  was  derived  for 
solution  by  digital  computer.  Thus  the  small  but  present  contribution 
from  distant  poles  and  zeros  were  taken  into  account. 

The  characteristic  determinant  of  the  stability  loop  was  formed  by 

the  coefficients  of  the  variables  defined  by  all  loop  transfer  functions. 

Referring  to  Fig.  2,  it  may  be  noted  that  these  variables  are  I 

N  or  r       Airframe  lateral  response 

N  ,  or  r.      Lateral  response  measured  by  instruments 

b_  Output  of  variable  gain  pot  and  compensation 

2  network,  if  included  1 

&d  Output  of  autopilot  actuators  % 

&o  Output  of  main  rudder  actuators 

Because  of  the  inability  of  the  IBM  704-  computer  to  accept,  in 

determinant  solution,  a  coefficient  of  higher  order  than  a  quadratic  the 

airframe  lateral  transfer  functions  were  not  used  as  such.  The  basic 

equations  ,of  motion,  for  three  degree  lateral  response,  from  which  the 

lateral  transfer  functions  were  derived,  were  substituted.  In  this  method 

all  coefficients  were  second  order  or  lower.  The  three  equations  of  motion 

are  developed  in  Chapter  II  of  Reference  10.   Assuming  rudder  deflections  as 

the  only  forcing  funrtion  and  after  application  of  perturbation  techniques 

and  the  LaPlace  transform  for  each  flight  condition  the  equations  take  the 

form: 
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(a^s*  ♦  a2s  +  a3)  p  ♦  (b,s2  ♦  b2s  ♦  b3)r  ♦  (c182  cas  c3)  <P    »  A  &u 

where  the  coefficients  of  the  variables  take  on  values  assigned  by 
aerodynamic  constants  and  derivatives  0  Any  of  the  constants  may  take 
the  value  of  zero  with  the  exception  of  the  constant  A. 

In  the  cases  where  N  is  the  lateral  response  being  plotted,  or 
some  portion  of  it  such  as  N '    *  r,  a  supplementary  equation  is  needed 
to  establish  the  relation  between  the  motion  equation  Vjariable  and  N  . 
This  equation  is  established  by  a  force  equation  in  the  lateral  direction 
which  results  in  J 

P       s    Y  ♦  mg  sin0    =    mV(p  ♦  r)  ♦  mVp 

0 

In  magnitude  comparison,  the  last  comparison,  the  last  term,  mVp,  is 
negligible  and  is  dropped.     Therefore,  the  lateral  force,  P   ,  becomes 
mV(p  ♦  r).     Then  it  follows  that  N    is  expressed  in  g's  by; 

N       s    !l     3    1   (P  ♦  *)  =  2   (*P  ♦  Tty 
7  mg  g  g 

where  it  conveniently  is  true  that  one  of  the  equations  of  motion  contains 
p  and  r  with  coefficients  of  1.0  and  a  substitution  to 'terms  of  p  and  fcu 
may  be  made. 

The  remaining  loop  variables  are  then  defined  and  interrelated  by  the 
transfer  functions  of  loop  components .  In  the  case  of  N  as  the  airframe 
response  and  filter  1  as  a  compensation  network,  these  equations  become  J 

Nyl(s*  +  175o9s  ♦  15795)  -  15795Ny  =  0 
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e 


Vo.is  - 1}  -  **  V  (B  * 1]  3  ° 

&o  (s2  ♦  40s  ♦  800)  -  fc^  (800)  =  0 

where  -  determines  feedback  sign  employed,  and  finally; 

^(s  ♦  10)  -  lO&o   =  0 

Therefore  the  coefficients  may  be  arranged  in  matrix  form  as  shown  in 
Fig.  A-l.  Thus  by  varying  the  value  of  K~_  and  solving  the  determinant, 
the  characteristic  equation  is  produced  and  may  be  factored  to  produce 
the  roots  coincident  with  the  particular  value  of  K-.  and  flight  condition 0 

In  a  similar  manner  the  matrices  for  the  other  responses  chosen  in  the 
analysis  can  be  derived .  The  number  of  loop  variables  present  then 
determining  the  number  of  simultaneous  equations  necessary  to  solve  the 
matrix  determinant.  The  control  equations  for  each  flight  condition  will 
not  vary  as  various  responses  are  tested.  However  the  response  N  does 
require  an  extra  equation  to  relate  this  force  to  the  control  variable 0 

In  summary  it  may  be  noted  that  the  characteristic  determinant  is 
actually  a  set  of  simultaneous  equations  whose  coefficients  determine  the 
characteristic  equation  of  the  loop  and  hence  its  roots.  It  is  emphasized 
again  that  due  only  to  computer  characteristics  were  the  basic  equations  of 
motion  substituted  for  the  airframe  transfer  functions  previously  given. 
The  control  equations  are  not  furnished  in  this  report  but  in  continuing 
work,  where  more  extensive  root  locus  investigation  is  necessary,  the 
transfer  functions  are  adequate  as  a  basis . 
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APPENDIX  B 
Analog  Computer  Implementation  and  Techniques 

As  a  part  of  this  study  analog  simulation  of  the  aircraft  and  control 
loop  was  made  and  computer  solutions  obtained.  This  study  was  carried  out 
for  two  basic  purposes.  First,  in  order  to  amplify  or  verify  the  root 
locus  characteristics  of  the  various  schemes  and  thereby  enable  the  choice 
of  the  most  acceptable  lateral  response.  Secondly,  in  order  to  provide 
the  capacity  for  proving  or  disproving  the  Adaptive  Controller's  capabilities. 
This  was  the  most  practical  way  to  test  the  adaptive  controller  in  the 
allotted  time. 

To  accomplish  the  first  part  of  the  analog  study  the  aircraft  control 
equations  and  equations  of  the  various  other  parts  of  the  system  were 
mechanized  on  Electronic  Associates  Analog  Computers.  This  mechanization  was 
done  in  accordance  with  the  diagram  in  Fig.  52.  The  lateral,  three  degree 
aircraft  control  equations  used  were  as  follows  J 

"t-   !l^  Cy  A   +  I*   MS^  -j>  +     h±    Cy.    SA 
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These  equations  of  motion  are  described  in  Chapter  II  of  Reference  10 
and  standard  symbols  are  employed.  In  the  mechanization  of  these  equations 
relays  were  used  to  permit  fast  switching  between  any  two  flight  conditions 
The  roll  channel  was  included  in  this  mechanization  in  the  hope  that  time 
would  permit  observation  of  the  effect  of  aileron  deflection,  6.  . 

The  transfer  functions  of  the  other  system  components  in  the  lateral 
loop  have  been  listed  in  the  text.  The  components  of  the  roll  loop  and 
their  transfer  functions  are  J 

Aileron  Actuators 

Roll  AFCS  Actuators 

Roll  Rate  Gyro        — r 

(s4  ♦  176s  ♦  15800) 

H2 


20 

-  55  deg 

(s  ♦  20) 
60 

(position) 
-  4-0  deg/sec 

(s  ♦  60) 
15800 

(rate) 

The  settings  of  the  potentiometers  for  the  various  flight  conditions 
are  tabulated  in  Table  XVIII.  The  flight  conditions  noted  on  these 
potentiometer  setting  sheets  correspond  to  conditions  1  through  6  as 
follows : 

Al Condition  1 

A2 Condition  2 

A3 Condition  3 

Bl Condition  A 

B2 Condition  5 

CI Condition  6 

Prior  to  data  recording,  the  steady  state  values  of  the  aircraft 
equations  were  checked.  Then  the  control  loops  were  closed  and  dynamic 
response  data  collected  at  each  flight  condition  for  each  response  analyzed. 

When  this  phase  of  the  analog  study  was  completed ,  the  Adaptive  Con- 
troller was  mechani  zed  in  order  to  test  its  performance  and  characteristics , 
The  diagram  of  this  mechanization,  as  well  as  a  switching  logic  table,  are 
contained  in  Fig.  53,  The  stepping  relay  timer  was  connected  to  and  oper- 
ated standard  eleven  pole  step  switches  from  a  digital  voltmeter.  This 
complete  mechanization  with  the  Adaptive  Controller  connected  to  modify 
the  variable  loop  gain  provided  the  facilities  to  completely  check  out  the 
logic,  operating  parameters,  and  aircraft  performance.  Recording  equipment 
was  included  as  needed  throughout  this  study. 
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